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ABSTRACT 


The  influence  of  a  solid  ice  cover  on  the  general  behaviour 
of  a  fixed  volume  oil  slick  retained  by  a  boom  has  been  studied.  A 
series  of  experiments  have  been  conducted  for  both  ice  free  and  ice 
covered  oil  slicks  encompassing  a  range  of  Froude  Numbers.  The  study 
addressed  deviations  in  profile  shape,  internal  oil  velocities,  and 
interfacial  waves  between  ice  free  and  ice  covered  oil  slicks.  An 
analytical  model  describing  the  internal  oil  velocity  distribution  has 
been  developed.  The  addition  of  an  ice  cover  was  found  to  increase 
the  oil  slick  thickness,  decrease  the  oil  slick  length,  and  suppress 
the  growth  of  interfacial  waves. 
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CHAPTER  1 


1.1  General 


INTRODUCTION 


During  the  past  decade,  considerable  investigative  research 
has  been  directed  towards  establishing  a  basic  understanding  of  the 
dynamic  behaviour  of  oil  slicks.  The  scope  of  these  previous  research 
studies  has  been  limited  to  the  occurrence  of  oil  spills  in  open  (ice 
free)  waters.  The  frequency  of  oil  spills  within  ice  covered  rivers 
will  probably  increase  as  the  development  of  northern  pipelines  is 
accelerated.  An  oil  spill,  in  all  likelihood  will  flow  under  an  ice 
cover  once  it  enters  a  river.  This  trend  dictates  the  necessity  for 
expanding  the  research  to  encompass  oil  spill  retention  under  ice 
covered  water.  The  objective  of  this  study  is  to  examine  particular 
characteri st ics  of  an  oil  slick  retained  by  a  boom  under  an  ice  cover 
in  flowing  water. 

1.2  General  Behaviour  of  Oil  Slicks  Retained  by  a  Boom 


Development  of  a  hypothesis  concerning  the  characteristics  of 

t 

oil  slicks  under  an  ice  cover  requires  prior  appreciation  of  the  basic 
behaviour  of  oil  slicks  in  open  flowing  waters.  For  this  study  the 
term  oil  slick  is  restricted  to  a  fixed  volume  of  oil  retained 
upstream  of  a  boom  within  a  fixed  width  of  channel.  The  following 
paragraphs  summarize  these  basic  characteri sties. 

Oil,  due  to  its  buoyancy,  floats  on  water  and  disperses  to 
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form  a  continuously  expanding  film.  The  spreading  of  an  oil  slick  is 
the  result  of  two  physical  forces.  First,  the  force  of  gravity 
induces  the  horizontal  flow  of  oil  due  to  the  excess  pressure  inside 
the  slick.  Second,  the  net  pull  of  surface  tension  forces  at  the  air 
-  water  -  oil  interface  acts  in  the  outward  radial  direction.  These 
two  driving  forces  are  opposed  by  inertial  and  viscous  forces.  This 
film  can  be  less  than  a  millimetre  thick. 

The  convective  movement  of  the  oil  slick  in  flowing  water 
coupled  with  the  extreme  thinness  of  the  oil  slick  dictates  that,  for 
economic  removal,  the  oil  must  be  ponded  to  a  sufficient  depth  by 
means  of  a  mechanical  barrier  termed  boom.  Any  attempt  to  draft  off 
the  oil,  without  an  increase  in  depth,  results  in  the  inclusion  of 
large  quantities  of  water,  thus  adding  to  the  time,  and  consequently, 
the  costs  of  removal. 

The  analysis  of  an  oil  slick  formed  by  towing  a  boom  in  calm 
water  is  similar  to  the  anlaysis  of  an  oil  slick  retained  by  a 
stationary  boom  in  flowing  waters,  since  the  ensuing  relative  current 
will  pond  the  oil  in  front  of  the  boom.  Oil  retained  by  a  boom 
assumes  a  distinctive  shape,  influenced  by  the  dynamic  and  viscous 
shear  forces,  as  a  result  of  the  relative  current  between  the  flowing 
water  and  arrested  oil.  Figure  1.1  depicts  the  characteristic  shape 
of  the  oil  slick  profile  at  three  different  phases  induced  by  various 
current  velocities. 

The  shape  of  the  arrested  oil  slick  has  been  divided  into 


3 


three  regions,  (indicated  in  Figure  1.1. b)  characterized  by  the 
relative  magnitudes  of  viscous  and  dynamic  forces  acting 
longitudinal ly. 

The  first  region  is  termed  the  headwave  or  frontal  zone  and 
portrays  the  shape  of  a  half  arrow-head.  The  headwave  is  a 
distinguishing  trait  of  all  two  phase  density  stratified  flows  in 
which  one  fluid  advances  on  the  other.  The  leading  edge  of  an  oil 
slick  is  similar  to  those  observed  when  salt  water  intrudes  into  fresh 
water,  when  a  cold  front  encroaches  on  a  warm  front,  and  when 
submerged  sediment  laden  water  flows  into  reservoirs.  The  shape  of 
the  headwave  is  governed  by  the  dynamic  force  exerted  on  the  oil  slick 
from  the  balance  between  horizontal  momentum  and  hydrostatic  force. 

The  second  region  is  designated  as  the  mid-region  or  viscous 
zone.  The  shape  of  the  viscous  zone  is  governed  by  the  cumulative 
viscous  shear  forces  which  supercede  the  dynamic  forces  after  a 
distance  of  about  10  -  20  thicknesses  from  the  leading  edge.  The 
viscous  shear  force  induces  a  "piling  up"  of  the  oil  which  when 
balanced  by  the  increasing  hydrostatic  pressure  in  the  oil,  results  in 
a  gradual  thickening  of  the  slick  from  its  minimum  thickness 
immediately  downstream  of  the  headwave. 

The  third  region  is  termed  the  near  boom  zone.  The  boom 
protrudes  deeper  than  the  oil  slick  forming  a  stagnation  point  as 
water  flows  beneath  the  boom,  which  induces  a  local  thinning  of  the 
oil  slick  upstream  of  the  boom. 
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The  prediction  of  equilibrium  oil  slick  thicknesses  for 
diverse  velocity  ranges  is  derived  from  detailed  analysis  of  the 
preceding  regions.  Typical  profiles  of  oil  slicks  at  three  such 
velocity  ranges  are  portrayed  in  Figure  1.  The  transition  between 
these  ranges  are  not  clearly  distinguishable  and  the  limits  vary  with 
the  properties  of  oils.  Order  of  magnitude  velocity  boundaries  are 
given  in  Figure  1.1  for  the  three  states  of  oil  slick  profiles. 

A  comprehensive  study  of  the  oil  slick  phenomenon  not  only 
considers  the  equilibrium  profile  described  above,  but  also  includes 
an  examination  of  oil  slick  behaviour  beyond  the  range  of  equilibrium 
velocity  limits.  Failure  of  a  boom  to  retain  oil  is  the  consequence 
of  two  independent  mechanisms.  The  first,  referred  to  as  drainage,  is 
simply  the  lack  of  sufficient  boom  draft  to  arrest  the  oil.  The 
thickness  of  the  slick  at  the  boom  is  predominantly  dictated  by  the 
relative  velocity  of  the  water  resulting  in  an  equilibrium  thickness 
and  length  of  slick  for  a  given  volume  of  oil  and  velocity. 

Therefore,  if  the  boom  draft  is  too  shallow,  the  oil  will  flow  under 
the  boom  sufficiently  reducing  its  volume  until  a  new  equilibrium 

thickness,  compatible  with  the  existing  boom  draft,  is  established. 

* 

Oil  slick  drainage  due  to  insufficient  boom  draft  occurs  very  quickly 
and  can  ultimately  result  in  considerable  loss  of  oil. 

The  second  mode  of  failure,  entitled  entrainment,  is 
independent  of  boom  design  and  occurs  at  higher  water  velocities. 
Superimposed  over  the  length  of  the  oil  slick  are  interfacial  waves 
travelling  longitudinally  from  the  leading  edge  of  the  oil  slick  to 
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the  boom.  When  two  different  liquids  move,  with  different  tangential 
velocities  across  their  interface,  the  kinetic  energy  of  motion  will 
cause  undulations  on  the  interface  to  grow  in  amplitude.  This 
phenomenon  is  called  Kelvin  -  Helmholtz  Instability. 

As  the  water  velocities  exceed  a  certain  critical  value, 
these  interfacial  waves  break,  resulting  in  oil  droplets  being 
entrained  into  the  flowing  water.  Owing  to  the  buoyant  properties  of 
the  oil,  these  droplets  will  rise  to  the  oil-water  interface  and 
migrate  towards  the  boom.  Consequently,  if  the  slick  is  of 
insufficient  length  so  that  the  oil  droplets  do  not  have  time  to  rise 
to  the  oil  slick  and  coalesce  before  reaching  the  near  boom  region, 
they  will  be  swept  under  the  boom.  It  is  the  loss  of  oil  which  does 
not  coalesce  and  is  swept  past  the  boom  that  constitutes  entrainment 
loss.  Even  though  the  critical  droplet  formation  velocity  is 
sustained,  there  may  be  no  actual  loss  of  oil.  The  majority  of  the 
wave  breaking  and  oil  droplet  entrainment  occurs  in  the  headwave 
region  as  indicated  in  Figure  1.1. c.  Prevention  of  oil  entrainment  is 
not  suppressed  by  a  well  adapted  boom  design  but  rather  by  the 
appropriate  deployment  of  the  boom. 

1.3  Scope  of  Study 

The  previously  mentioned  description  of  the  characteri stic 
traits  of  oil  slick  behaviour  in  open  waters  are  well  documented  in 
available  literature  (Wicks  (1969),  Lindenmuth,  Miller,  and  Hsu 
(1970),  Cross  and  Hoult  (1971),  Wilkinson  (1971,  1973),  Hale,  Norton, 
and  Rodenberger  (1974),  and  Milgram  and  Van  Houten  (1978)).  However, 
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the  current  state  of  knowledge  regarding  the  effects  of  an  ice  cover 
is  deficient.  This  thesis  studies  the  differences  in  behaviour  of  ice 
covered  oil  slicks  from  ice  free  oil  slicks. 

The  scope  is  restricted  to  the  treatment  of  a  solid  ice  cover 
as  a  boundary  condition  and  does  not  consider  deviations  in  dynamic 
behaviour  as  induced  by  changes  in  the  physical  properties  of  the  oil 
resulting  from  interaction  between  oil  and  ice.  The  emphasis  of  the 
study  addresses  the  basic  characteristics  of  the  oil  slick  profile  and 
interfacial  stability  rather  than  the  detailed  design  of  booms  or 
recovery  methods.  As  indicated  earlier,  the  oil  slick  under  question 
is  composed  of  a  fixed  volume  of  oil  retained  upstream  of  a  rigid  boom 
in  a  fixed  width  flowing  channel. 

The  approach  used  in  the  remainder  of  this  study  is  as 
follows.  First,  a  presentation  of  a  detailed  literature  review  of  ice 
free  oil  slicks  is  given  in  Chapter  2.  Chapter  3  outlines  the 
laboratory  study  of  oil  slicks  both  with  and  without  an  ice  cover. 
Chapter  4  develops  various  mathematical  relationships  describing  the 
expected  differences  between  ice  free  and  ice  covered  oil  slicks. 
Chapter  5  presents  the  experimental  results  and  Chapter  6  develops  an 
interpretation  of  these  results  which  leads  to  an  explanation  for  the 
increased  thickness  of  the  oil  slick  profile  under  an  ice  cover  and 
for  the  apparent  stability  of  the  oil-water  interface.  Finally, 
Chapter  7  summarizes  the  conclusions  and  suggests  recommendation  for 
further  study. 
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CHAPTER  2 


LITERATURE  REVIEW 


2.1  General 

The  objective  of  this  study,  (as  previously  stated),  is  to 
examine  specific  characteristics  of  an  oil  slick  retained  by  a  boom 
under  an  ice  cover  in  flowing  water.  As  a  preparatory  phase,  an 
extensive  survey  of  existing  literature  was  conducted.  This  survey 
encompassed  the  analytical  and  experimental  studies  pertaining  to  the 
headwave  phenomena,  the  effects  of  currents  on  viscous  zone 
development,  oil  entrainment  and  drainage  failure,  interfacial 
stability  and  other  related  subjects.  The  major  emphasis  of  the 
reviewed  studies  was  restricted,  in  scope,  to  ice  free  conditions, 
whereas  the  limited  amount  of  literature  addressing  the  problem  of  ice 
covers  neglected  the  dynamic  behaviour  of  oil,  but  rather  investigated 
the  physical  properties  of  the  oil  and  their  interaction  with  ice. 

This  chapter  presents  a  concise  summary  of  the 
afore-mentioned  oil  slick  studies  through  the  examination  of 
individual  components  involved  in  the  oil  slick  phenomena. 

2.2  Headwave 


The  leading  edge  of  all  density  stratified  flows  is 
customarily  designated  as  the  headwave  (Figure  2.1. a).  The 
characteristic  shape,  and  subsequently  the  analysis,  of  an  oil  slick 
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headwave  is  analogous  to  that  of  submerged  gravity  currents 
investigated  by  several  researchers  as  saltwater  intrusions. 
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Dimensional  analysis  reveals  that  the  details  pertaining  to 
stratified  flow  can  be  described  by  three  basic  non-dimensional 
parameters.  These  parameters  are:  the  Reynolds  Number  (ratio  of 
inertial  to  viscous  forces),  the  Weber  Number  (ratio  of  inertial  to 
interfacial  tension  forces),  and  the  Densimetric  Froude  Number  (ratio 
of  inertial  to  buoyant  forces).  Wilkinson  (1971)  has  demonstrated 
that  flow  in  the  vicinity  of  the  headwave  is  dominated  by  the  dynamic 
and  buoyant  forces.  Consequently  the  effects  of  interfacial  surface 
tension  and  viscosity  are  secondary,  thus  the  headwave  phenomena  can 
be  adequately  characterized  by  the  Densimetric  Froude  Number. 
Designating  the  approach  velocity  as  U,  the  specific  gravity  of  oil  as 
s,  with  the  headwave  thickness  as  h,  then  the  Densimetric  Froude 
Number  (hereafter  termed  the  Froude  Number)  is  defined  as: 

F  (2.2.1) 

[gd-sjh]1'2 

Von  Karman  (1940) 

* 

An  ideal  fluid  model  for  steady  propagation  of  gravity 
currents  derived  from  the  application  of  the  steady  state  irrotational 
Bernoulli  Theorem,  applied  between  the  stagnation  point  forward  of  the 
leading  edge  of  the  headwave,  and  points  on  the  interface  downstream 
was  postulated  by  Von  Karman.  Assuming  that  the  interface  becomes 
horizontal,  the  resulting  relationship  between  the  propagation 
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velocity  (U)  of  the  gravity  current  and  current  thickness  (h)  is: 


h  =  U2  (2.2.2) 

2g  ( 1— s ) 

In  comparison  to  Equation  (2.2.1)  this  yields  a  Froude  Number  of  yj 2  . 

Benjamin  (1968) 


Benjamin  presents  a  comprehensive  treatment  involving  the 
application  of  inviscid  fluid  theory  to  steady  gravity  currents,  in 
which  he  disputes  Von  Karman's  assumption  of  energy  conservation  as 
unjustifiable.  Rather,  Benjamin  states  that  a  more  fundamental 
condition  is  the  overall  balance  of  momentum  flux  against  hydrostatic 
forces  in  the  fluid. 

The  theoretical  analysis  of  the  headwave,  regardless  of 
whether  the  gravity  current  is  the  bottom  fluid  (as  in  the  case  of  oil 
slicks)  is  similar.  In  this  analysis,  the  gravitational  constant  g 
will  appear  multiplied  by  a  similar  factor  of  (1-s).  Expanding  this 
concept,  Benjamin  argues  that  the  upper  fluid  could  be  replaced  by  an 

f 

air-filled  cavity  displacing  liquid  beneath  a  horizontal  boundary. 
However,  to  establish  continuity  within  this  review,  the  following 
theoretical  results  extacted  from  Benjamin's  paper  include  the  factor 
(1-s). 


Equating  the  change  in  the  momentum  flux  with  the  pressure 
forces  between  an  approaching  section  of  depth  D  and  a  downstream 
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section  beneath  the  oil  slick  of  depth  d,  coupled  with  the  equation  of 
continuity,  the  following  relationship  for  the  downsteam  depth  d  as  a 
function  of  the  approach  velocity  U  is  obtained. 


U2  =  d  (D2  -  d2)  (2.2.3) 

g  ( 1 — s )  0  D2  (2D  -  d) 

Alternately,  expressed  in  terms  of  the  headwave  thickness  h  =  D-d, 
Equation  (2.2.3)  defines  the  Froude  Number  as: 


F  =  (D  -  h)  (2D  -  h)  (2.2.4) 

h  D  (D  +  h) 

This  definition  of  the  Froude  Number  for  the  headwave  is 
derived  on  the  basis  of  momentum  balance. 

A  second  equation  incorporating  a  headloss  h^  can  be 
obtained  from  the  application  of  the  Bernoulli  energy  equation  between 
the  same  approaching  section  and  receding  downstream  section.  The 
momentum  equation  combined  with  this  energy  equation  yields  the 
following  relationship  for  energy  loss  hL  in  terms  of  the  headwave 
thickness: 


(D  -  2h)  h2 

2(D2  -  h2) 


(2.2.5) 


If  energy  was  conserved,  (i.e.  hL=0),  as  assumed  by  V on 
Karman,  then  the  non-trivial  solution  of  Equation  (2.2.5)  reveals  that 
the  headwave  thickness  h  must  occupy  exactly  one  half  the  depth  of 
flow  in  the  approaching  section.  However,  the  headwave  thickness,  as 
reported  by  various  investigators,  assumes  a  range  of  values  all  of 
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which  are  considerably  less  that  one  half  the  approaching  flow  depth. 
Moreover,  observed  headwaves  are  characterized  by  the  leading  front 
possessing  a  wave  like  form  (so  called  "head")  whose  peak  obtains  a 
thickness  up  to  twice  the  downstream  height.  Observations  reveal  that 
the  leeward  side  of  this  wave  is  a  highly  turbulent  zone  suggesting 
some  form  of  wave  breaking.  This  wave  breaking  process  is  indicative 
of  an  energy  loss  mechanism,  and  Equation  (2.2.5)  dictates  the  energy 
loss  mechanism,  and  Equation  (2.2.5)  dictates  the  energy  loss 
requirements,  for  various  headwave  thicknesses,  necessary  to  satisfy 
the  governing  momentum  Equation  (2.2.4). 

Specifying  hm  as  the  maximum  thickness  of  the  leading  wave 
trough,  and  using  the  concepts  of  irrotational  flow  beneath  the 
receding  oil  slick  (at  thickness  h)  is  subcritical.  This  condition  is 
satisfied  when  the  range  of  downstream  thicknesses  exceeds 
h  /  D  >  0.6527. 

The  headwave  observed  at  the  front  of  oil  slicks  is 
essential,  not  merely  incidental,  to  the  dissipative  processes 
necessary  to  satisfy  both  the  momentum  and  energy  considerations. 
Utilizing  an  integral  form  of  momentum  balance  and  Bernoulli's 
theorem,  Benjamin  established  that  magnitude  of  the  maximum  headwave 
thickness  hm  must  adhere  to  the  following  relationship: 

hm  /  h  >  (2.2.6) 

For  the  limit  of  infinite  depth.  Equation  (2.2.4)  yields  the  value  of 
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the  Froude  Number  F^  as  yj 2  .  Keulegan  (1958)  reports  on  his 

experiments  involving  saltwater  wedges  a  ratio  of  hm/h  as  2.36  which 
complies  with  the  theoretical  condition  of  Equation  (2.2.6). 

Wilkinson  (1971) 


The  previous  theoretical  analysis  and  accompanying 
experimental  studies  were  motivated  by  the  need  to  investigate  gravity 
currents  other  than  the  oil  slick  phenomena.  A  direct  analysis  of  the 
oil  slick  headwave  (which  closely  parallels  the  arguments  employed  by 
Benjamin)  was  conducted  by  Wilkinson.  As  previously  disclosed, 
Wilkinson  demonstrated  the  dominance  of  dynamic  over  viscous  forces  in 
the  vicinity  of  the  headwave.  The  momentum  flux  under  the  oil  slick 
results  in  the  slick  experiencing  a  balancing  pressure  force  of  the 

O 

order  p  Ui.  The  slick  is  simultaneously  subjected  to  a  viscous 
shear  force,  resulting  from  the  cumulative  interfacial  shear  stress  of 
the  order: 


The  ratio  of  shear  forces  to  dynamic  pressure  forces  at  any  section 
from  the  front  of  the  oil  slick  will  be  of  the  order: 

2  h 

Extracting  from  the  experimental  work  of  Cross  and  Hoult  (1971),  the 
value  of  the  interfacial  shear  coefficient  of  an  oil  slick  C.  as 
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0.01,  Wilkinson  concludes  that  the  dynamic  forces  exceed  the  viscous 
forces  by  a  factor  of  10  for  approximatly  20  slick  thicknesses 
downstream  from  the  front  of  the  slick. 

Equating  the  dynamic  pressure  forces  with  the  momentum  flux, 
Wilkinson's  resulting  momentum  equation  for  the  Froude  Number  in  terms 
of  a  relative  depth  H  =  h/D  is: 

Fh  =  [  |  2 H  +  1,_1  ]  1/2  (2.2.7) 

'1-H  s ' 

Examining  the  case  of  infinite  depth.  Equation  (2.2.7)  yields 
the  value  of  the  Froude  Number  F^  as  (2s)  The  resulting  value 
for  the  headwave  thickness  h  deviated  from  Benjamin's  prediction  by 
the  factor  1/s.  This  difference  stems  from  the  replacement  of 
Benjamin's  air  cavity  under  a  solid  boundary  with  an  oil  slick  at  a 
free  surface.  Wilkinson's  analysis  accounts  for  the  rise  of  the  oil 
slick  above  the  mean  free  surface  of  the  approaching  section.  If 
Wilkinson's  stagnation  point  was  taken  at  the  top  of  the  oil  surface, 
then  Benjamin's  results  are  obtained. 

Experiments  on  peanut  oil,  conducted  by  Wilkinson, 
corroborate  the  theoretical  relationship  of  Froude  Number  variation 
with  relative  depth.  Comparison  between  various  other  experimental 
studies  and  the  theoretical  analysis  presented  by  Benjamin  and 
Wilkinson  reveals  a  tentative  verification  of  the  ideal  fluid  theory 
but  also  exhibits  deviations  arising  from  real  fluid  properties. 
Keulegan  reported  a  Froude  number  of  F^  =  1.2  for  a  saltwater  wedge 


' 
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with  a  relative  depth  ratio  of  h/D  of  0.111.  The  corresponding 
theoretical  Froude  Number  of  -  1.23,  as  predicted  from  Equation 
(2.1.4)  is  in  close  agreement.  Oil  slick  studies  performed  by 
Lindenmuth,  Miller,  and  Hsu  (1970)  describe  values  of  the  Froude 

Number  F^  ranging  from  a  low  of  1.0  to  a  high  of  1.3.  Lindenmuth  et 
al  noted  that  the  Froude  Number  decreased  as  the  oil  slick  thickness 
(and  consequently  the  relative  depth)  increased,  which  is  a  direct 
interpretation  of  the  afore-mentioned  theoretical  analysis.  However, 
consideration  of  the  experimental  facilities  utilized  indicates  that 
there  was  insufficient  variation  in  the  relative  depth  to  account  for 
the  range  of  observed  Froude  Numbers.  Two  real  fluid  properties  were 
considered  to  be  instrumental  in  producing  the  variations  in  Froude 
Numbers.  It  was  remarked  that  increasing  interfacial  surface  tension 
tends  to  reduce  the  thickness  of  the  headwave  at  a  given  speed. 

A  second  observed  property  is  that  an  increase  in  viscosity 
apparently  induces  a  greater  thickness  of  the  headwave  at  a  given 
water  current. 

Benjamin  resolved  in  his  discussion  of  experiments  by  Zukoski 
(1966)  involving  the  advance  of  an  air  bubble  in  a  long  water  filled 
tube,  that  the  reduction  in  observed  propagation  velocity,  as  noted  by 
Zukoski,  were  smaller  than  predicted  by  theory  due  to  the  effects  of 
surface  tension.  The  surface  tension  influenced  the  propagation 
velocity  by  two  distinct  mechanisms.  First,  the  neglect  of  surface 
tension  would  over  estimate  the  balancing  pressure  in  the  headwave  at 
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the  stagnation  point,  and  secondly,  the  downstream  forces  balancing 
the  momentum  flux  should  be  lower  due  to  the  contribution  of  surface 
tension  at  the  interface.  Though  Zukoski's  experiments  revealed  that 
surface  tension  effects  resulted  in  a  reduction  of  the  propagation 
velocity,  an  alternative  interpretation,  similar  to  that  of  the 
Lindenmuth  et  al  study,  is  that  the  surface  tension  contributes  to  a 
reduction  in  the  headwave  thickness  at  a  given  water  current. 

Hale,  Norton,  Rodenberger  (1974) 


Realizing  the  role  of  real  fluid  properties,  Hale  et  al, 
developed  an  equation  for  the  headwave  thickness  in  which  the  effects 
of  surface  tension  (T)  and  viscosity  (y  )  were  i ncorporated. 
Correlating  the  Froude  number  with  Reynolds  and  Weber  numbers,  by  the 
utilization  of  a  least  squares  fit  to  their  experimental  data,  the 
following  relationship  for  the  maximum  headwave  thickness  (h  )  was 
developed ; 


,  7  „  3.254  0.0157 

hm  =  7. 535. U  «w 

7,  ,1.325  t0.325 

g  (l-s)  t 


(2.2.8) 


This  relationship  indicates  the  magnitude  of  the  influence 
the  effects  of  real  fluid  properties  exercise  on  the  development  of 
the  headwave  thickness.  Similar  to  the  observations  of  Lindenmuth  et 
al,  Equation  (2.1.8)  reveals  that  comparatively,  as  the  viscosity 
increases,  so  does  the  headwave  thickness,  and  that  increases  in 
surface  tension  result  in  a  decrease  in  headwave  thickness. 
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Milgram  and  Van  Houten  (1978) 


In  a  recent  study,  Milgram  and  Van  Houten  dispute  the 
hypothesis,  as  prescribed  by  Wilkinson  and  others,  that  the  dynamic 
forces  predominate  the  viscous  forces  in  determining  the  headwave 
prof i le. 


The  afore-mentioned  analysis  by  Benjamin  and,  again,  by 
Wilkinson,  indicated  that  an  energy  loss  must  occur  to  satisfy  the 
momentum  flux  requirements.  If  viscous  shear  forces  are  neglected  the 
only  mechanism  for  dissipation  of  energy  is  the  breaking  of  the 
leeward  side  of  the  headwave  which  results  in  the  generation  of 
rotational  flow  and  a  consumption  of  energy.  However,  Milgram  and  Van 
Houten  advocate  that  the  consideration  of  a  free  surface  and  shear 
stress  at  the  oil  water  interface  remove  the  theoretical  necessity  for 
breaking  of  the  headwave. 

Various  experimental  studies  including  Wicks  (1968), 

Wilkinson  and  Lindenmuth  et  al  all  observed  that  the  headwave  would 
break  only  at  the  higher  velocities  tested.  Their  studies  focused 
solely  on  the  spatially  constant  shear  stress  in  the  viscous 
downstream  zone,  and  as  a  result  neglected  to  include  the  effects  of 
shear  stress  in  the  vicinity  of  the  headwave. 

The  major  emphasis  of  the  Milgram  and  Van  Houten 
investigation  was  the  disclosure  of  the  relative  importance  of  shear 
stress  not  only  in  the  downstream  viscous  zone  but,  also,  its  role  in 
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the  development  of  the  headwave  profile.  Evaluation  of  the  oil  slick 
profile  and  its  dependence  on  shear  stress  was  accomplished  through 
the  following  analysis. 

It  was  assumed  that  a  potential  flow  approximation  of  the 
flowing  water  phase  was  accurate  for  predicting  pressures  in  the 
water.  Thus  the  pressure  along  the  bottom  of  the  oil  layer  was 
calculated  for  potential  flow  along  the  oil  slick  profile  measured 
experimentally.  This  calculated  pressure  distribution  and  the 
measured  oil  thickness  distribution  were  in  turn  implemented  to 
determine  the  shear  stress  distribution  utilizing  an  equation  relating 
the  interfacial  shape,  pressure  and  shear  stress.  Although  the 
complete  details  of  the  development  of  the  equation  are  contained  in  a 
following  section  dealing  with  the  viscous  zone,  the  basic  rationale 
of  the  equation  consists  of  a  force  balance  for  a  differential  element 
of  the  oil  slick.  In  order  to  compare  the  relative  effects  of  friction 
and  pressure  on  the  determination  of  the  interfacial  shape,  a  new  oil 
slick  thickness  distribution  in  the  absence  of  dynamic  forces  was 
calculated  by  the  previous  equation  using  the  calculated  shear 
stresses  and  setting  the  dynamic  pressure  and  its  horizontal 
derivatives  to  zero.  Milgram  and  Van  Houten  present  several  plots  of 
the  experimentally  measured  profiles  and  calculated  profiles  based  on 
shear  stress  for  a  complete  range  of  experimental  conditions.  In  the 
range  of  lower  velocities,  and  consequently  for  oil  slicks  without  the 
observed  breaking  headwave,  the  two  profiles  match  satisfactorily. 

The  most  remarkable  feature  of  the  analysis  was  that  the  calculated 
friction  coefficients  in  the  vicintity  of  the  headwave  assumed  a  range 
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of  values  exceeding  those  previously  reported  in  the  literature  by  a 
factor  of  10.  As  a  result  of  the  apparent  matching  of  the  calculated 
and  measured  oil  slick  profiles  and  large  shear  at  the  headwave 
Milgram  and  Van  Houten  concluded  that  the  previously  cited  assumption 
that  friction  is  unimportant  in  the  headwave  development  is 
ill-founded.  However,  observation  of  the  calculated  profiles  at  the 
high  velocities  revealed  that  their  analysis  deteriorates  to  the  point 
of  inaccurately  describing  the  development  of  a  breaking  headwave. 

2.3  Viscous  Zone 


In  the  mid  region  of  the  oil  slick  between  the  headwave  and 
the  restraining  boom  the  oil  slick  undergoes  more  gradual  thickening 
than  in  the  headwave.  Contrary  to  the  headwave  zone  where  the  dynamic 
pressure  forces  contribute  to  the  development  of  the  oil  slick 
profile,  the  viscous  shear  stresses  are  predomi nantly  responsible  for 
the  profile  in  the  mid  region.  Hence,  this  region  is  entitled  the 
viscous  zone  (Figure  2. 1. b) . 

The  viscous  shear  forces  are  cumulative  along  the  length  of 
the  slick  downstream  from  the  leading  edge  of  the  oil,  thus 
contributing  to  a  gradual  increase  in  the  thickness  of  the  oil  slick. 
The  cumulative  shear  forces  are  balanced  by  an  increasing  hydrostatic 
pressure  gradient  corresponding  to  the  increasing  thickness  of  oil. 

Associated  with  the  viscous  shear  forces  are  interfacial 
waves  superimposed  on  the  mean  oil  slick  profile.  These  waves 
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originate  at  the  leading  edge  of  the  oil  slick  and  travel  downstream 
through  the  viscous  zone  with  a  characteristic  shape  and  velocity. 
Section  2.4  (Interfacial  Instabilities)  covers  this  phenomenon  in 
detail.  The  remainder  of  this  section  is  concerned  strictly  with  the 
development  of  the  mean  oil  slick  profile  in  the  viscous  zone. 

Cross  and  Hoult  (1971) 

One  of  the  original  analysis  of  the  viscous  zone  was 
presented  by  Cross  and  Hoult.  Their  experiments  were  conducted  for  a 
low  range  of  water  velocities  at  which  no  noticeable  headwave  was 
formed  and  thus  their  analysis  did  not  distinguish  between  a  headwave 
and  a  viscous  zone. 

Assuming  the  momentum  flux  within  the  oil  to  be  an  order  of 
magnitude  less  than  the  momentum  in  water.  Cross  and  Hoult  equated  the 
hydrostatic  pressure  gradient  arising  fron  the  depth  of  oil  h,  to  the 
interfacial  shear  stress  This  force  balance  takes  the  form  of; 

p  g(l-s)  h  dh  =  a,  =  1  p  C,  U2  (2.3.1) 

"dx  1  7  w  1 

If  the  coefficient  of  friction  C^  is  assumed  to  be 
constant,  as  suggested  by  their  experiments,  then  the  intergration  of 
Equation  (2.3.1)  yields  the  following  relationship  for  the  oil  slick 
prof i le; 


1^1 

. 
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l>2  =  u2  C,  x  (2.3.2) 

gTT^)  1 

The  oil  slick  profile  assumes  a  parabolic  shape  with  the 
thickness  growing  to  the  one  half  power  of  the  longitudinal  distance. 
The  experimental  profiles  revealed  values  for  the  coefficient  of 
friction  to  be  0.005  for  No.  2  oil  and  0.008  for  soybean  oil. 

Cross  and  Hoult's  argument  for  a  spacial  constant  value  of 
the  coefficient  of  shear  friction  C.#  used  to  evalute  the 
interfacial  shear  as  a  function  of  the  approach  water  velocity,  were 
based  on  the  following  analogy  employing  a  sand  roughened  flat  plate. 
At  low  Reynolds  Number  Rx  (Reynolds  Number  expressed  with 
longitudinal  distance  x  as  a  length  scale),  the  plate  appears  smooth 
and  C.  will  decrease  with  increasing  distance  x.  However,  at  higher 
values  of  Rx  (supposedly  after  the  headwave  region)  if  the  grain 
roughness  K$  increases  wih  distance  such  that  the  ratio  Ks/x 
remains  constant  then  win  attain  a  constant  value.  Cross  and 
Hoult  observed  that  the  interfacial  waves  would  grow  in  magnitude  as 
they  migrated  downstream,  and  thus,  they  argued  that  the  increasing 
wave  amplitudes  provided  an  increasing  roughness  height  Ks, 
sufficient  to  yield  a  constant  value  for  the  coefficient  of  friction. 

Wicks  (1969) 


A  considerably  more  complex  analysis  of  the  viscous  zone  was 
developed  by  Wicks  in  which  he  utilized  experiments  encompassing  a 
broader  range  of  flow  conditions  conducted  on  a  much  larger  scale. 
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Wicks  observed  the  headwave  zone  to  be  distinct  from  the  viscous 
zone.  Based  on  a  concept  of  boundary  layer  separation  and 
reattachment,  Wicks  defined  the  start  of  the  viscous  zone  at  a 
distance  of  25  hm  (maximum  headwave  thickness)  from  the  leading  edge 
of  the  oil  slick.  Wicks  draws  an  analogy  to  Plate's  (1964) 
investigation  of  abrupt  surface  roughness  in  a  turbulent  boundary 
layer,  in  which  reattachment  occurred  about  50  disturbance  heights 
downstream.  Benjamin  has  shown  that  the  maximum  headwave  thickness  to 
be  approximately  twice  the  thickness  of  the  downstream  neck  of  the 
headwave.  Wicks  took  the  difference  in  thickness  to  be  analogous  to 
an  abrupt  surface  roughness  and  consequently  assumed  that  reattachment 
occurs  at  a  distance  of  25  times  the  maximum  headwave  thickness  as 
given  by  Benjamin. 

Reattachment  of  the  boundary  layer  in  flowing  water  induces  a 
drag  on  the  oil  interface  which  in  turn  sets  up  a  circulation  within 
the  oil  layer.  Wicks'  analysis  of  the  progressive  thickening  of  the 
oil  slick  involves  the  simultaneous  solution  of  the  following  three 
equations.  First,  the  circulation  in  the  oil  phase  is  related  to  the 
interfacial  drag  by  the  application  of  the  Navier-Stokes  equation. 
Assuming  a  linear  velocity  distribution  between  the  downstream 
oil-water  interfacial  velocity  and  the  upstream  oil-air 
interfacial  velocity  U^.,  will  yield  a  relationship  between  the 
oil-water  interfacial  velocity  Ui  and  the  interfacial  shear  stress 
a.j  as  fol  lows; 


' 
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(2.3.3) 


Boundary  layer  theory  is  used  to  describe  the  local  shear 


stress  caused  by  the  water  flow  as  a  function  of  the  difference 
between  the  free  stream  velocity  and  interfacial  velocity  as; 


°1  -  1/2  P„  C-j  (U  -  U^2 


(2.3.4) 


Wicks  adopts,  from  Schlichting,  two  relationships  for  the 


coefficient  of  friction  as  functions  of  the  Reynolds  Number  R 

A 

Finally,  a  force  balance  in  the  horizontal  direction  shows  that  the 
slope  of  the  interface  must  exactly  balance  the  interfacial  shear 
stress ; 

a.  =  pn  g  h  dh  (2.3.5) 

1  0  dx 

Numerical  solution  of  these  relationships  yields  a  prediction 
of  the  slick  thickness  as  a  function  of  the  longitudinal  distance. 


A  literature  review  by  Lau  and  Kirchefer  (1974)  outlines  two 


criticisms  of  Wicks'  analysis.  First,  the  assumption  of  a  linear 
velocity  profile  within  the  oil  is  inconsistent  with  other  reported 
observations  in  which  the  oil-water  interfacial  velocity  exceeds  the 
oil  surface  velocity  several  times.  Secondly,  Wicks'  force  balance 
Equation  (2.3.5)  has  neglected  the  component  of  the  hydrostatic 
pressure  in  the  direction  of  the  sloping  interface. 
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Wilkinson  (1971,  1973) 


A  third  method  of  analysis  for  the  viscous  zone  was  described 
by  Wilkinson.  The  emphasis  of  his  study  focused  on  the  endeavor  to 
incorporate  the  boundary  effects  inherent  in  a  finite  depth  of  flow. 
Wilkinson's  force  analysis  of  the  headwave  zone  (Section  2.1) 
indicated  that  after  some  20  thicknesses  downstream  any  further 
increase  in  thickness  was  due  to  viscous  effects.  Downstream  of  the 
headwave  zone,  Wilkinson  utilized  the  following  three  assumptions  to 
develop  a  viscous  zone  model: 

1)  The  momentum  flux  due  to  circulation  in  the  oil  slick  is 
small  in  comparison  to  the  interfacial  shear  forces. 

2)  The  flow  beneath  the  oil  slick  is  steady  and  uniform  in 
distribution  over  the  depth  of  flow. 

3)  The  equilibrium  of  the  retained  oil  slick  dictates  that  any 
change  in  the  longitudinal  gradient  of  inertial  and  pressure 
forces  must  be  balanced  by  the  boundary  shear  stress. 

This  third  premise  equally  applies  to  a  control  element 
composed  of  the  oil  slick  and  underlaying  water  flow,  or  a  control 
element  of  the  oil  slick  alone  (Figure  2.2). 

Wilkinson  derived  the  following  equation  for  an  overall 
momentum  balance  in  terms  of  unit  discharge  q; 

d  |-low  g  s  (d+h)2+  1  PW  g(l-s)  d2  +  pwq2  ]  =  -ob 
dx  ’  3 


(2.3.6) 


■ 
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Consideration  of  the  forces  acting  on  a  control  element 
consisting  of  the  oil  slick  alone  yields; 

d  j-l_p0  g  s  h2]  =  „.  -p0  g  s  h  1(d)  (2.3.7) 

dx  2  dx 

The  boundary  and  interfacial  shear  stresses  are  defined 
respectively  as; 

ab  =  1/2  C5  pw  U2  and  ^  =  1/2  Ci  pw  U2 

In  terms  of  the  densimetric  Froude  Number,  Equations  (2.3.6) 
and  (2.3.7)  can  be  combined  non-dimensional ly  to  form; 

2H  dH  [D*  -  (H  +  ,  F  ^)]  =  ,  F  ^  ["  Term  A  1  (2.3.8) 

dX  V|VH;  Dt~H; 

Term  A  =  C,  [D.-(l-s)  (H  -  ,J_  ^2)]  +  C  H 

—  D. -H'  D 

s  t 

In  which  H  =  h  ,  D.  =d  ,  X  =  x  and  F  =  U 

D  D  “  [g ( 1-s )D] 1/2 

Equation  (2.3.8)  contains  a  singularity  for  the  interfacial 
slope  dH/dX  when; 


Dt  =  H  + 


<Dt-H> 


(2.3.9) 


At  lesser  values  of  H,  the  interfacial  slope  dH/dX  is 
positive,  indicating  a  thickening  of  the  oil  slick.  The  interfacial 
slope  becomes  infinte  for  the  singularity  at  the  critical  thickness. 
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Finally,  for  larger  values  of  H  the  interfacial  slope  is  negative. 

A  peculiar  chacteristic  of  the  critical  oil  slick  thickness, 
given  by  Equation  (2.3.9),  is  its  independence  from  viscous  forces. 
However,  the  maximum  length  of  a  slick  at  which  this  critical 
thickness  occurs  depends  totally  on  the  magnitude  of  viscous  stresses. 


Numerical  integration  of  Equation  (2.3.9),  between  the  oil 
slick  thickness  at  the  neck  of  the  headwave  and  the  critical 
thickness,  indicates  the  maximum  length  and  consequently  the  volume  of 
the  viscous  zone  is  a  function  of  the  Froude  Number  and  the  magnitude 
of  the  interfacial  (C^)  and  boundary  (C^)  coefficients  of 
friction.  This  procedure  reveals  that  the  dimensionless  depth  Dj 
does  not  exceed  unity  by  more  than  a  few  percent  and,  therefore,  the 
maximun  dimensionless  slick  thickness  can  be  adequately  expressed  as; 


H  =  1-F2/3 


(2.3.10) 


Small  scale  experiments  performed  by  Wilkinson  indicate  the 
interfacial  shear  stress  coefficient  C.  is  not  spacially  constant  as 
originally  assumed,  but  rather,  decreases  with  the  longitudinal 
distance.  This  concept  was  observed  as  a  more  gradual  thickening  of 
the  oil  slick  than  that  predicted  by  Equation  (2.3.8).  If  Wilkinson 
had  incorporated  a  spacially  varying  coefficient  of  friction  into  his 
analysis,  the  singularity  of  Equation  (2.3.8)  and  resulting  critical 
thickness  would  not  have  been  affected,  since  the  critical  thickness 
criteria  stems  from  a  Froude  Number  instability. 
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Lindenmuth,  Miller  and  Hsu  (1970) 


In  an  attempt  to  investigate  the  role  of  waves  on  the 
performance  of  various  oil  slick  retention  devices,  a  series  of  large 
scale  oil  slick  tests  were  conducted  by  Lindenmuth  et  al.  A 
prelimenary  set  of  experiments  were  performed  in  the  absence  of  waves 
to  establish  a  set  of  control  data. 

In  an  effort  to  support  the  experimental  findings,  a 
simplified  analysis  of  the  viscous  zone  was  derived  by  the 
consideration  of  a  differential  control  volume  of  oil.  The  net 
hydrostatic  forces  acting  on  the  control  were  equated  with  the  shear 
stress  at  the  interface  to  yield  the  following 

h  dh  =  sum  (2.3.11) 

dx  pq  g(l-s) 

This  force  balance  of  the  control  volume  differs  from  that 
given  by  Cross  and  Hoult  in  Equation  (2.3.1)  by  the  factor  pQ  /  p 
The  difference  is  due  to  the  allowance  for  the  oil  slick  to  float 

above  the  water  free  surface.  The  slick  thickness  h  is  related  to  the 

* 

depth  of  oil  below  the  water  surface  m  by; 


Lindenmuth  et  al  realized  that  the  major  difficulty  in 
applying  this  theory  lay  in  the  complexity  of  evaluating  the  shear 


- 


. 
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forces.  Available  boundary  layer  theories  are  inadequate  to  account 
for  the  effects  of  interfacial  waves,  flow  separation  at  the  headwave, 
oil  viscosity,  and  circulation  within  the  oil  layer.  However,  it  is 
suggested  that  for  engineering  estimates,  adequate  oil  slick 
geometries  could  be  calculated  on  the  assumption  of  a  constant  shear 
force  along  the  slick.  Integration  of  Equation  (2.3.11)  coupled  with 
a  constant  shear  assumption  yields  a  parabolic  shape  of  the  oil  slick 
thickness  distribution  in  the  following  form; 

h(x)2  =  G ( x )  +  h^  (2.3.12) 

2 

in  which  G  =  U  and  hn=  thickness  at  the  neck  of  the 
pQ  g(l-s)  headwave. 

Application  of  Equation  (2.3.12)  to  the  experimental  data 
reveals  a  range  of  interfacial  shear  coefficients  from  0.002  for 
diesel  fuel  to  0.013  for  freshwater  over  saltwater.  Engineering 
estimates  of  the  volume  of  oil  retained  by  a  boom  is  given  by  the 
integration  of  Equation  (2.3.12); 

V(x)  =  Vh  +  2/3G  (hfj  +  G(x)  )2/3  -  (2.3.13) 

Based  on  experimental  data  for  the  headwave  geometries  the 
volume  of  oil  in  the  headwave  is  given  by; 


Vh  =  15.6  hfj  (ft2/ft) 


■ 
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Lindenmuth  et  al  propsed  that  utilization  of  Equations 
(2.3.12)  and  (2.3.13)  would  provide  sufficient  estimation  for  the 
required  boom  draft  to  prevent  drainage  failure  of  a  given  quantity  of 
oil,  in  the  absence  of  waves. 

Although  this  study  is  not  directly  concerned  with  the 
effects  of  approaching  waves  on  oil  slicks,  the  following  comments  on 
the  development  of  the  viscous  zone,  influenced  by  waves,  as  noted  by 
Lindenmuth  et  al,  are  briefly  discussed.  The  purpose  of  their 
presentation  here  is  to  provide  a  comparative  background  for  the  brief 
discussion  of  the  detailed  study  of  wave  interaction  on  interfacial 
velocities  and  shear  stresses  outlined  in  the  next  section. 

The  range  of  tests  conducted  indicated  that  under  certain 
conditons  a  significant  increase  in  the  oil  slick  thickness  could 
occur  due  to  waves.  Typically  the  slick  would  thicken  at  the  wave 
crests  and  decrease  in  the  waves  troughs.  The  development  of  the 
viscous  zone  appeared  to  be  a  complex  function  of  the  wave  heights, 
wave  period,  and  boom  response  to  waves.  Lindenmuth  et  al  were  unable 
to  provide  a  quantitative  theory  as  to  their  relative  roles  but 
several  trends  were  noted.  The  wave  steepness,  defined  as  the  ratio 
of  wave  height  to  wave  length,  was  found  to  be  the  predominant 
parameter  influencing  the  oil  slick  thickness.  The  ratio  of  mean 
slick  thickness  in  waves  to  mean  slick  thickness  in  calm  water  was 
found  to  be  non-linearly  proportional  to  wave  steepness.  The  maximum 
increase  in  slick  thickness  occurred  when  the  incident  waves  were 
breaking. 
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The  headwave  was  strongly  disturbed  under  certain  conditions, 
suggesting  that  the  phenomena  does  not  result  from  a  simple 
superimposition  of  velocities.  Furthermore, it  was  noted  that  the 
large  increase  in  slick  thickness  occurred  for  the  full  range  of 
frequencies  tested,  indicating  that  a  simple  resonance  phenomenon  was 
not  the  responsible  mechanism.  It  was  discerned  that  when  the  wave 
frequency  was  sufficiently  high,  such  that  wave  lenghts  were  less  than 
3  times  the  boom  draft,  these  waves  were  completely  reflected  by  the 
boom  with  no  resulting  increase  in  slick  thickness. 

In  conclusion,  Lindenmuth  et  al  suggest  that  due  to  the 
complex  nature  of  waves  effects  on  the  viscous  zone,  any  boom  draft 
requirements  dictated  by  a  "current  only"  analysis  may  have  to  be 
increased  by  a  factor  of  2  or  more. 

Hale,  Norton,  and  Rodenberger  (1974) 


A  most  comprehensive  investigation  of  the  effents  of  currents 
and  waves  on  an  oil  slick  was  conducted  by  Hale,  Norton  and 
Rodenberger.  The  major  emphasis  of  their  study  was  directed  towards 

3 

the  development  of  a  mathematical  oil  droplet  entrainment  model,  and 
towards  the  evaluation  of  various  oil  slick  recovery  devices. 

However,  extensive  preliminary  study  of  the  interface  within  the 
viscous  zone  was  conducted  to  establish  a  complete  understanding  of 
the  oil  slick  phenomenon.  This  study  consisted  primarily  of  the 
collection  of  a  comprehensive  set  of  velocity  data;  the  resulting 
implications  are  presented  briefly  in  the  following  section.  Details 


; 


. 
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concerning  the  precise  magnitude  for  the  data  are  presented  later  in 
this  study  to  complement  velocity  data  acquired  by  author  during  the 
course  of  the  experimental  phase  of  the  investigation. 


Using  a  pitot  tube,  at  several  different  sections  along  the 
oil  slick,  velocity  profiles  in  the  water  below  the  oil  slick  were 
established.  A  non-linear  least  squares  regression  was  employed  to 
match  the  observed  velocities  to  the  following  distribution; 

U  =  _  A2  exp  “^3^  (2.3.14) 

Where  y  is  the  distance  below  the  free  surface  and  A^s  ^  and  A3 
are  the  computed  coefficients. 

This  form  of  velocity  distribution  was  in  agreement  with  the 
observed  velocities.  Although  the  velocity  data  was  applied  to  a 
one-seventh  power  law  velocity  distribution,  it  was  unable  to 
duplicate  the  same  degree  of  correlation  as  in  the  previous 
distribution. 

Defining  the  shear  stress  at  the  interface  as; 


Hale  et  al  evaluated  the  shear  stress  by  the  differentiating 
Equation  (2.3.14).  This  process  lead  to  an  interesting  observation 
concerning  the  magnitude  of  shear  stress  at  various  sections  along  the 
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slick.  Hale  et  al  reported  that  the  interfacial  shear  coefficients  in 
the  vicinity  of  the  headwave  were  1/100  of  those  in  the  mid-region  as 
reported  by  Cross  and  Hoult.  This  led  Hale  et  al  to  the  conclusion 
that  Wilkinson's  assumption  of  the  insignificance  of  shear  stresses  in 
the  frontal  zone  is  valid. 

Expanding  on  the  concept  of  Cross  and  Hoult  that  the 
interfacial  shear  stress  is  balanced  by  the  horizontal  hydrostatic 
pressure  gradient  corresponding  to  the  increase  in  the  thickness  of 
oil  (Equation  (2.3.1)),  Hale  et  al  developed  the  following 
relationship  between  the  average  interfacial  shear  coedfficient  C-, 
the  average  thickness  hav  and  the  slick  length  1; 

C.  =  9  g(l-s)  h2  (2.3.15) 

i  a  '  av  v  ' 

4  1  U2 

Substitution  of  the  experimental  data  into  Equation  (2.3.15) 
reveals  relatively  the  same  shear  coefficients  as  previously  reported 
by  Cross  and  Hoult. 

In  addition  to  the  water  velocity  measurements,  an  attempt 
was  made  to  develop  a  model  for  the  velocities  within  the  oil  slick 
itself.  The  velocity  gradient  in  the  oil  at  the  interface  was  derived 
by  equating  the  shear  stress  in  the  water  to  that  in  the  oil.  Thus, 
the  oil  velocity  gradient  at  the  interface  is; 
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dU 


i  oil 


y  dU 
—  —  ) 

yo  ^  i  water 


An  additional  boundary  condition  was  established  by  timing 
the  velocities  of  chips  on  the  oil-air  surface.  These  two 
measurements  coupled  with  the  known  zero  shear  on  the  surface  and  the 
requirement  of  mass  continuity  within  the  circulating  oil  layer  lead 
to  the  development  of  a  velocity  profile  in  the  oil  layer.  The 
coefficients  of  the  following  relationship  were  evaluated  for  these 
considerations ; 


U  =  B 


1 


+  BoYB3_1  exp  -B4yB3 


(2.3.16) 


Where  Y  is  the  dimensionless  depth  (ie.  0.0  at  the  free 
surface,  and  1.0  at  the  interface) 


However,  due  to  the  random  nature  of  interfacial  velocities 
and  gradients  arising  from  the  complication  introduced  by  interfacial 
waves,  the  results  of  this  analysis  were  considered  to  be  questionable 
by  Hale  et  al. 

* 

A  dimensionless  interfacial  velocity  (U^/u  )  relationship 

was  developed  by  equating  the  velocities  to  the  Reynold  Number  Re> 

viscosity  ratio  (y  /  y  )  and  the  dimensionless  slick  position 

u  w 

(x/t).  Using  a  least  square  procedure,  the  following  equation  was 
deri ved ; 


■ 
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U.  =  6.20 


-0.0397 


-0.133 


R 


-0.214 


(2.3.17) 


e 


Though  the  scatter  in  data  was  considerable,  Equation 
(2.3.17)  indicates  a  strong  trend  for  a  decrease  in  U^,  with 
increasing  x  (the  downstream  distance). 

A  major  emphasis  of  the  investigation  by  Hale  et  al  was 
directed  towards  the  study  of  the  combined  effects  of  currents  and 
waves.  Although  the  details  of  their  study  are  beyond  the  scope  of 
this  study,  some  of  the  basic  conclusions  regarding  the  interfacial 
shear  stresses  and  velocities,  as  affected  by  waves,  have  been 
extracted  and  are  presented  briefly  in  the  following  paragraphs. 


The  oil  slick  set-up  induced  by  waves  is  mainly  an  effect  of 


the  mass  transport  in  waves  and  energy  transfer  from  waves  to  the  oil 
slick.  The  most  simple  analysis  is  accomplished  by  the 
superimposition  of  a  progressing  wave  system  velocity  field  on  the 
velocities  of  the  current  alone.  Hale  et  al  have  demonstrated  that  if 
the  phase  shift  between  the  motions  in  oil  and  water  is  zero,  then  the 
interfacial  shear  will  be  independent  of  the  waves.  However,  if  a 
phase  shift  occurs,  the  shear  stress  takes  the  following  form; 


=  1/2  Ci  Pw(U2  +  q2w2(l-costf ) ) 


(2.3.18) 


a  ■ 


Their  analysis  is  insufficient  to  predict  the  magnitude  of 
phase  shift,  but  is  does  provide  a  mechanism  by  which  waves  can  alter 
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the  interfacial  shear  and  consequently  the  oil  slick  length. 

A  model  test  utilizing  a  compliant  barrier  (wave  refeclts 
less  than  10  percent  of  incident  wave)  were  conducted  over  a  complete 
range  of  velocities  and  wave  climates.  At  low  current  velocities, 
waves  tended  to  increase  the  slick  length  with  the  exception  of  steep 
waves.  At  higher  current  velocities  all  waves  induced  a  decrease  of 
the  oil  slick.  The  periodic  water  velocities  caused  local  variation 
in  slick  thickness  in  a  manner  similar  to  that  observed  by  Lindemuth 
et  al.  Hale  et  al  were  unable  to  describe  the  effects  of  waves  in  any 
systematic  manner.  The  only  conclusion  postulated  is  that  waves  could 
significantly  alter  slick  lengths  but  their  simple  analysis  of  the 
interfacial  shear  was  inadequate  to  fully  describe  this  complex 
process. 

Milgram  and  Van  Houten  (1978) 

An  ingenious  analysis  of  the  oil  slick  problem  is  presented 
in  a  recent  paper  by  Milgram  and  Van  Houten.  The  prime  objective  of 
their  analysis  was  to  evalute  the  relative  importance  of  interfacial 
shear  stress  and  dynamic  pressure  in  determining  the  thickness 
distribution  of  an  oil  slick.  As  mentioned  in  the  previous  section 
dealing  with  headwaves,  Milgram  and  Van  Houten  argue  that  if  shear 
stresses  in  the  frontal  zone  are  sufficiently  high  then  the 
theoretical  necessity  for  a  breaking  headwave,  as  suggested  by 
Benjamin,  is  removed.  Milgram  and  Van  Houten  cite  as  preliminary 
evidence,  that  the  concept  that  dynamic  forces  are  not  dominant  in  the 
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headwave  region,  one  of  the  basic  observations  derived  during 
experimental  work  with  oil  slicks.  More  explicitly,  the  leading  edge 
of  the  slick,  as  observed  by  all  authors,  is  an  order  of  magnitide 
less  than  60  degrees  to  the  horizontal,  as  predicted  by  the  dynamic 
analysis  of  Von  Karman. 

Their  analysis  and  subsequent  arguments  supporting  the 
postulation  of  equal  dominance  of  shear  stresses  and  dynamic  forces 
stems  from  the  following  combination  of  theory  and  experiment.  In 
theory,  it  was  assumed  that  potential  flow  approximates  the  flow 
within  the  water  phase  to  a  sufficient  degree  of  accuracy.  Thus, 
utilizing  experimental  measurements  of  the  oil  slick  thickness  as  a 
boundary  condition,  numerical  analysis  could  calculate  the  flow  and 
corresponding  pressure  distribution.  A  theoretical  model  equating  the 
pressures,  oil  thickness  shape,  and  shear  stresses  was  then  utilized 
to  evaluate  the  shear  stress  distribution  from  the  measured  oil  shape 
and  calculated  pressure  distribution. 

Numerical  integration  of  the  water  velocity  potential 
provided  the  tangential  velocity  Vt  along  the  measured  oil  slick 
shape.  The  dynamic  pressure  P  along  the  interface  was  then  given  by 
the  Bernoulli  equation,  in  a  non-dimensional  form  as; 

p  =  1  [1  -  (  Ut  )2] 

U 


(2.3.19) 


. 
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Arguing  that  measurements  by  Hale  et  al  indicate  that  the 
velocity  in  the  oil  are  sufficiently  smaller  than  those  in  the  water, 
Milgram  and  Van  Houten  assumed  that  dynamic  pressure  and  momentum 
fluxes  within  the  oil  can  be  neglected.  The  remaining  mean  forces  are 
shown  in  Figure  2.3. a.  These  are,  pressure  forces  on  the  vertical 
faces,  stready  and  dynamic  pressure  forces  along  the  lower  interface, 
and  the  combined  viscous  shear  and  Reynolds  stresses  along  the 
interface.  The  interfacial  stresses  can  be  expressed  in  terms  of  the 
approach  velocity  U  and  a  spacial  variant  coefficient  of  friction 

C.j.  Pressure  continuity  relates  the  mean  thickness  of  oil  h  to  the 
oil  thickness  h  as  measured  below  the  original  free  surface  in  the 
following  form; 

h  =  pw  g  m  +  P  (2.3.20) 

^~~g 

The  measured  slick  profile  m  and  the  horizontal  length  scale 
x  can  be  described  non-dimensional ly  by  the  factor  g/U^. 

Such  that  M  =  m  g  and  L  =  x  g 

U2  U2 

A  balance  of  the  horizontal  forces  on  the  differential 
element  of  oil,  coupled  with  the  afore-mentioned  factors,  yield  the 
following  basic  equation  relating  the  shape  of  the  oil  to  the  dynamic 
pressure  and  shear  forces  at  the  interface. 

(M  +  P)[(l-s)dM  +  dP]  =  1  s  C.  (2.3.21) 

dL  dL  2  1 
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The  experimental  portion  of  Milgram  and  Van  Houten's  study 
encompassed  a  broad  range  of  water  velocities  and  oil  properties. 

They  have  presented  several  plots  of  these  various  experiments  which 
consisted  of  the  measured  interface  shape  m,  the  calculated  pressure 
distribution  P  and  the  calculated  shear  stress  distribution  C-.  In 
order  to  evaluate  the  relative  effects  of  friction  and  pressure  in 
determining  the  interfacial  shape,  Milgarm  and  Van  Houten  developed  a 
relationship  to  predict  the  oil  thickness  distribution  in  the  absence 
of  dynamic  pressures.  This  relationship  consisted  of  setting  the 
dynamic  pressure  and  its  horizontal  derivative  to  zero  in  Equation 
(2.3.21)  and  solving  for  M.  This  solution  was  defined  as  M^; 

M^(x)  =  (  s  C - ( x )  dx)  ^  (2.3.22) 

1-s 

This  frictional  oil  thickness  distribution  was  plotted  with 
the  above  mentioned  distributions.  Figure  2.3. b  and  2.3.C  show  two 
typical  plots  as  traced  from  Milgram  and  Van  Houten.  These  two  plots 
reveal  the  two  distinct  shapes  of  the  calculated  interfacial  friction 
distributions.  The  Type  I  distribution  (Figure  2. 3. b)  was  observed  to 
have  a  maximum  shear  in  the  vicinity  of  maximum  headwave  thickness. 

The  magnitude  of  the  maximum  friction  coefficient  ranged  from  0.013  to 
0.073  which  is  an  order  of  magnitude  greater  than  the  friction 
coefficients  previously  reported  by  other  investigators.  The  friction 
distribution  decreases  to  a  range  between  0.004  to  0.012  behind  the 
headwave.  Type  I  distributions  were  observed  when  the  interfacial 
waves  were  diminutive.  However,  when  the  interfacial  waves  exceeded  a 
specific  value  the  shear  distribution  took  another  shape  as  indicated 
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in  Figure  2.3.C.  The  Type  II  distribution  were  characterized  by  more 
than  one  peak  in  the  maximum  value  of  friction  coefficients.  The 
first  peak  occurred  upstream  of  the  maximum  oil  thickness  in  the 
headwave  with  a  second  peak  immediately  downstream  of  the  headwave. 
The  magnitude  of  the  coefficients  for  Type  II  distributions  ranged  up 
to  0.33. 


The  calculated  oil  thickness  was  observed  to  positively 
match  the  measured  oil  thickness  M  for  Type  I  distributions.  The  two 
plots  were  practically  identical  in  the  viscous  zone.  Based  on  their 
calculations  Milgram  and  Van  Houten  state  the  assumption  that  the 
friction  was  instrumental  in  determining  the  geometry  of  the  rear 
portion  of  the  oil  layer  propounded  by  several  previous  investigators, 
is  justified.  However,  due  to  the  apparent  match  of  the  two  oil 
thicknesses  in  the  vicinity  of  the  headwave,  Milgram  and  Van  Houten 
conclude  that  the  previous  assumption,  specifically  that  the  fricion 
is  unimportant  for  the  forward  portion  of  the  layer,  is  ill-founded. 

2.4  Interfacial  Instability 


The  preceding  two  sections  have  dealt  with  the  equilibrium 
characteristics  of  an  arrested  oil  slick.  However,  the  successful 
containment  of  an  oil  slick  must  include  consideration  of  the  various 
aspects  of  boom  failure.  Aside  from  the  obvious  failure  resulting 
from  insufficient  boom  draft,  termed  oil  drainage,  a  complex  failure 
mechanism  entitled  entrainment  failure  must,  also,  be  taken  into 
consideration.  Entrainment  failure  stems  from  the  interfacial  shear 


, 
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instabilities  and  the  resultant  formation  of  interfacial  waves.  In 
the  vicinity  of  the  headwave,  the  interfacial  waves  break  and  oil 
droplets  are  entrained  in  the  water.  The  actual  magnitude  of  oil  loss 
by  this  mechanism  is  dependent  on  the  ability  of  the  oil  droplets, 
once  entrained  into  the  flowing  water,  to  rise  back  to  the  slick 
interface  and  coalesce  before  being  swept  under  the  boom.  This 
entrainment  failure,  which  occurs  at  current  velocities  of 
approximately  0.4  meters  per  second,  is  invariant  to  the  actual  design 
of  the  boom  and,  thus  becomes  a  fundamentally  limiting  factor  in  oil 
spill  retention  and  recovery. 

The  development  of  numerous  models  for  predicting  the  oil 
droplet  entrainment  losses  have  resulted  from  the  various  studies  of 
oil  slick  entrainment  failure.  These  models,  though  different  in 
degree  of  complexity,  consist  primarily  of  an  evaluation  of  the  rates 
of  droplet  formation  at  velocities  exceeding  the  critical  formation 
velocity,  examination  of  the  droplet  charateristics,  trajectories  and 
terminal  velocities,  and  an  investigation  of  the  reattachment 
phenomena.  The  details  of  these  entrainment  loss  models  is  beyond  the 
scope  of  the  experimental  portion  of  this  study.  Thus,  the  remainder 
of  this  section  will  briefly  examine  various  falure  criteria  and  then 
addressed  the  development  of  interfacial  waves  resulting  from  K-H 
instability.  This  treatment  shall  consist  of  an  examination  of  the 
basic  problem  and  method  of  analysis  of  K-H  instability  waves  with 
emphasis  on  its  application  to  oil  slick  interfaces. 


. 
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Wicks  (1969) 


Since  the  oil  droplet  entrainment  phenomena  was  first 
reported  by  Wicks  (1969),  several  theoretical  and  empirical 
correlations  for  the  critical  velocity  Ucr,  of  entrainment 
initiation,  and  various  models  of  entrainment  losses  have  been  evolved. 

Wicks  suggested  that  droplet  entrainment  is  inevitable  as  a 
result  of  the  unstable  headwave  predicted  by  Benjamin's  analysis. 

Wicks  observed  interfacial  waves  formed  at  the  leading  edge  of  the  oil 
slick,  which  grew  in  magnitude  as  they  migrated  downstream  toward  the 
headwave.  Droplets  were  torn  from  the  rear  of  the  headwave  and 
entrained  in  the  flow.  Based  on  experimental  observations,  Wicks 
stated  that  the  velocity  at  which  the  oil  droplets  are  torn  from  the 
headwave  appears  to  be  in  agreement  with  the  predictions  of  a  critical 
Weber  Number,  as  suggested  by  Hinze,  for  gradual  acceleration  of  the 
dispersed  phase.  In  terms  of  the  critical  velocity  Ucr,  the  average 
droplet  diameter  d,  and  the  interfacial  surface  tension  T,  the  Weber 
Number  is  defined  as; 

(We>  =  0W  ucr  d  =  14-°  (2.4.1) 

T 

Wicks  further  stated  that  the  formation  of  waves  at  the 
oil-water  interface,  within  the  viscous  zone,  corresponded  well  with 
the  predictions  of  Kelvin  -  Helmholtz  instability  phenomena  (hereafter 
termed  K-H  instability). 
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Lindenmuth,  Miller,  and  Hsu  (1970) 


Lindenmuth  et  al  proposed  that  droplet  formation  behind  the 
headwave  is  due  to  form  drag.  As  water  flows  under  the  headwave,  a 
low  pressure  region  develops  aft  of  the  headwave.  Then,  as  the 
velocity  of  the  current  increases,  the  low  pressure  eventually  becomes 
sufficient  to  tear  oil  away  from  the  rear  of  the  headwave,  with  the 
resulting  formation  of  oil  droplets.  Experimentally,  the  critical 
velocity  at  which  oil  droplet  formation  commences  was  observed  to  be 
primarily  a  function  of  interfacial  surface  tension. 

Wilkinson  (1973) 


Wilkinson  indicated  that  the  development  of  interfacial  waves 
was  in  some  repects  similiar  to  the  growth  of  wind  induced  waves  on  a 
free  surface.  The  occurence  of  interfacial  waves  suggested  that 
energy  was  being  transferred  from  the  main  flow  into  the  oil  slick, 
possibly  statisfying  the  energy  dissipation  requirements  of  Benjamin's 
headwave  theory.  However,  Wilkinson's  experimental  observations  of 
the  interfacial  wave  characteristics  revealed  that  only  five  percent 
of  the  available  energy  was  consumed  by  the  flux  of  wave  energy. 

Wilkinson  had  realized  that  the  K-H  instability  was 
responsible  for  the  production  of  interfacial  waves,  but  he  made  no 
reference  to  its  role  as  a  mechanism  for  the  entrainment  phenomena. 
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Hale,  Norton,  and  Rodenberger  (1974) 

The  extensive  investigation  conducted  by  Hale  et  al  suggested 
that  the  use  of  a  critical  Weber  Number  is  an  adequate  means  to 
evaluate  the  critical  velocity  for  oil  droplet  formation.  Their 
experimental  work  confirmed  that  droplet  formation  occurs  strictly  at 
the  aft  of  the  headwave  but  that  actual  entrainment  of  the  oil 
droplets  exists  throughout  the  entire  length  of  the  slick. 

Kelvin  -  Helmholtz  Instability 


The  basic  principle  behind  investigations  of  the  stability  of 
fluid  flow  is  that  given  an  initial  flow,  the  superimposition  of  small 
perturbations  upon  the  flow  are  examined  to  determine  whether  the 
amplitude  of  the  disturbances  increase  or  decrease  with  time,  thus 
determining  whether  a  given  flow  is  stable  or  unstable. 

K-H  type  of  instability  is  that  which  occurs  when  two 
different  liquids  move  with  different  tangential  velocities  across 
their  interface.  The  momentum  exchange  induces  any  wave  disturbance 
on  the  interface  to  grow  in  amplitude  to  the  verge  of  mixing  the  two 
liquids.  Interfacial  surface  tension  tends  to  reduce  the  unstable 
growth,  as  does  viscosity,  but  to  a  secondary  degree. 

Alternately,  K-H  instability  can  be  interpreted  as  the 
resulting  action  of  the  pressure  distribution  in  phase  with  the 
interface  elevation,  attempting  to  overcome  the  "stiffness"  (due  to 
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the  restoring  force  of  gravity  and  surface  tension)  of  the  interface. 
This  type  of  interpretation  is  typically  followed  by  an  examination  of 
the  instability  of  a  vortex  sheet  at  the  oil-water  interface. 

Though  rigid  mathematical  analysis  of  vortex  sheets  at  liquid 
interfaces  has  been  studied  in  some  detail,  it  will  be  the  former 
description  of  the  K-H  instability  phenomena  that  will  be  addressed  in 
the  remaining  paragraphs. 

The  basic  mathematical  analysis  of  K-H  instability  is  well 
documented  in  existing  literature.  The  succeeding  analysis  follows 
closely  that  presented  by  Lamb,  (Art  232). 

A  simple  case  of  instability  occurrs  when  two  fluids, 
homogenous  within  themselves,  of  densities  pw  and  Pq,  one  beneath 
the  other  (pw  >  Po),  move  in  a  parallel  direction  with  uniform 
velocities  Uw  and  U0.  The  definition  sketch  of  Figure  2.4  is 
developed  for  three  phase  flow  (addressed  later  in  this  section)  but, 
also,  serves  as  a  definition  sketch  for  two  phase  flow.  If  the  flow 

in  the  upper  and  lower  layers  is  assumed  to  be  irrotational ,  then  the 

* 

velocity  potentials  of  <bQ  and  dw  which  satisfy  the  Lapace  Equation 
exist. 


A2  d0  =  0  and  A2  =  0  (2.4.2) 

If  the  velocity  potentials  are  disturbed  by  infinitesimally 
small  perturbation  dQ'  and  tfw'  so  that  ; 
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\  =  Uox  +  K  and  K  =  Uwx  +  K'  (2.4.3) 

It  follows  that  the  perturbations  themselves  must,  also, 
adhere  to  the  Laplace  equations. 

a2  =  0  and  a2  tfw'  =  0  (2.4.4) 

Considering  the  interface  between  the  two  liquids  to  be  the  y 
origin  and  assuming  the  liquids  infinite  in  both  directions,  then  the 
dynamic  boundary  condition  (B.C.)  for  the  perturbation  is; 

lim  y  »  infinity  =  0  and  (2.4.5) 

1 im  y  >  -infinity  d  '  =  0 

Utilizing  these  two  B.C.s  the  solution  to  Equation  (2.4.4) 
can  take  the  following  form; 

$0‘  =  AQ  exp(-Ky  +  i(Kx  -  at))  (2.4.6) 

V  =  Aw  exP(  +  i(Kx  "  at))  (2.4.7) 

In  which  a  =  KC  and  K  =  wave  number  (positive  and  real) 

C  =  wave  velocity  (complex) 

The  interface  between  the  two  liquids  is  a  streamline  of  both 
fluids  and  the  path  of  a  particle  on  the  interface  is  described  as 
y  =  n(x, t )  . 
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The  total  derivative 

D_  (y  -  n(x, t) )  =  0  (2.4.8) 

Dt 

is  a  kinematic  boundary  condition  which  indicates  that  once  a  particle 
is  on  the  interface,  it  stays  on  the  interface.  Expanding  the  total 
derivative  (Equation  (2.4.8))  and  realizing  that  vertical  velocity  V 
can  be  expressed  in  terms  of  the  velocity  potential  results  in  two 
additional  kinematic  B.C.'s. 


For  the  top  fluid; 


dn  +  U  dn  =  -d01  at  y  =  n 

dt  0  ?t  -r-5. 

dy 

and  for  the  bottom  fluid; 


dn.  +  U  dn  =  -d0 '  at  y  =  n 
dt  w  dt  — 

dy 


(2.4.9) 


(2.4.10) 


If  the  amplitude  of  n  is  assumed  to  be  small,  then 
utilization  of  Taylor's  Expansion  allows  the  previous  two  B.C.  to  be 
evaluated  at  y  =  0. 


An  additional  dynamic  boundary  condition  results  from 

* 

consideration  of  the  continuity  of  pressure  at  the  interface.  (In 
this  case  surface  tension  is  neglected.)  Recalling  that  the  interface 
is  a  streamline,  the  Bernoulli  Equation  for  the  top  fluid  is; 


The  symbol  d  represents  the  partial  derivative 
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d0'  +  1  (U^  +  vj:  )  +  P  +  gn  =  C(constant) 
3T  7  0  0  o 


or  to  a  first  order  approximation. 


Po  =  -  Uo  ~  +  C 


(2.4.11) 


"dT 


W~ 


Similarily  the  Bernoulli  Equation  for  the  bottom  fluid  is; 


P  =  -d0 '  -  U  d0'  -  gn  +  C 
w  w  ai 


(2.4.12) 


w 


3x 


The  two  fluids  share  a  common  streamline  at  the  interface  and 
the  pressure  P  is  constant  on  both  sides  of  this  streamline  (in  the 
absence  of  surface  tension.)  Thus,  equating  equation  (2.4.11)  and 
(2.4.12)  yields; 

po  (d0o  +  Uo  %  +  9n)  =  pw  (d0w  +  Uw  +  9n)  (2.4.13) 

dt  dx  dt  dx 

Defining  the  form  of  the  interface  n  as; 


n  =  a  exp[i (Kx  -at)] 


(2.4.14) 


Utilizing  the  general  solutions  for  the  perturbations  0  ■ 

and  0w«  (ie.  Equations  (2.4.6)  and  (2.4.7)),  the  dynamic  boundary 
condition  based  on  continuity  of  pressure  at  the  interface  (equation 

(2.4.13))  becomes; 

°o(iAo  p  +  Uo’KAo  +  ga)  =  °w(_iV  +  UwiKAw  +  ga)  (2.4.15) 

The  kinematic  boundary  conditions  at  the  interface  (Equation 
(2.4.9)  and  (2.4.10))  can  be  reduced  in  a  similar  manner  by 
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substituting  Equation  (2.4.6),  (2.4.7),  and  (2.4.14)  to  form; 


For  the  top  fluid; 


ia(U0K  -  „)  =  A0K 


(2.4.16) 


and  for  the  bottom  fluid; 


1a(UwK  -  a)  =  -V 


(2.4.17) 


The  simultaneous  solution  of  Equations  (2.4.15),  (2.4.16), 
and  (2.4.17)  for  the  three  unknows  AQ,  Aw,  and  a  yield; 

e0(KU0-a)2  +  Pw(KD„-a)2  =gK(pw-p0)  (2.4.18) 


The 

velocity  C  = 


quadratic  solution  of  Equation  (2.4.18)  for  the  phase 
o/K  ,  of  an  interfacial  disturbance  is; 


a 

K 


p  U  +p  U  i  [g  p  -p 
Mo  o  w  w  Mw  Mo 

p  +  p  K  p  +p 

o  w  w  o 


P  P  (U  -U  )2]1/2 
(p  +p  )2 


(2.4.19) 


The  first  term  is  referred  to  as  the  weighted  (by  density) 
mean  velocity  of  the  two  fluid  streams.  Relative  to  this  velocity, 
interfacial  waves  travel  with  a  wave  velocity  ±C  as  given  by; 


-  p  p 

Mo  w 


(p  +p  ) 


(Uo- 


w 


2 


(2.4.20 


. 
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2 

Where  CA  =  a  p  -  p 
0  y 


CQ  is  the  wave  velocity  of  disturbances  in  the  absence  of 
currents.  These  waves  would  subside  with  time  and  result  in  a  stable 

flow.  However,  the  interfacial  waves  would  be  unstable  if  C  is 
imaginary.  That  is,  if; 

<Uo  -  Uw>2  >  9  >w  -  Pp  (2-4.21) 

K  po  pw 

The  results  of  this  form  of  analysis  reveal  that  in  the 
absence  of  any  stabilizing  forces,  such  as  surface  tension,  the 
interface  will  be  unstable  for  any  small  disturbance. 

The  neglect  of  interfacial  surface  tension  between  oil  and 
water,  is  totally  unjustified  for  the  evaluation  of  K-H  instabilities 
of  oil  slicks.  Incorporating  the  effects  of  surface  tension  into  the 
above  analysis  will  only  alter  the  dynamic  boundary  condition  of 
pressure  continuity  at  the  interface.  The  previous  analysis  assumed 
that  pressure  at  the  interface  of  the  two  fluids  were  equal.  However, 
the  inclusion  of  surface  tension  dictates  that  the  surface  tension  a  , 
must  be  balanced  by  the  pressure  difference  at  the  interface:  ie. 


-T  (d2ri  +  d2n) 
dx2  dy2 


at  y  =  o 


(2.4.22) 


Applying  the  two  linearized  Bernoulli  Equations  (2.4.11)  and 
(2.4.12)  to  Equation  (2.4.22)  at  y  =  o  yields; 
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po(d^  +W  +9n)  -°w(d#lw  +Uwd,lw  +9n)  ■  -T  d2n  (2.4.23) 
dt  dx  dt  dx  d? 

Pursuing  the  similiar  solution  of  the  previous  example,  the 
resulting  second  order  polynomial  for  wave  velocity  C  and  the 
following  roots; 


C  =  pwUw  +  poUo  *  tc02  -  P0PW  (Up  -  UW)2  ]1/2 


0  +  0 
MW  H0 


+  P 


W 


In  which  C ^  is  now  C  =  g  pul  -  p^  +  T  K 


w 


K  pw  +  po 


D  +  0 
MW  M0 


(2.4.24) 

(2.4.25) 


As  previously  indicated  the  first  term  of  Equation  (2.4.24) 
is  the  weighted  mean  velocity  and  CQ  is  the  interfacial  velocity  of 
disturbance  waves  in  the  abence  of  currents.  Thus  as  before  if  C  is 
positive  and  imaginary,  instabi 1 ity  wi 11  occur,  ie.  if; 


p 


o 


w 


>  C 


2 

0  min 


(2.4.26) 


Cq  is  a  function  of  the  disturbance  wave  number  K  ,  and  has 
a  minimum  value  of; 


C0  min  “  [  2Tg  (pw  -  p„)  ]  1/2  (2-4.27) 

(po  +  pw)2 

at  a  critical  wave  number  K  of; 

Kcr2  =  9  (pw  -  p0)  (2-4.28) 


T 
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Therefore,  surface  tension  will  suppress  the  K-H  instability 
below  a  critical  velocity  difference  of; 

<u0  -  uw>2  <  2t9  (pw  -Oq)  (2.4.29) 

po  pw 

#■ 
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Jones  (1972) 


The  previous  K-H  instability  analysis  deals  with  two 
semi-infinite  inviscid  fluids  of  different  density,  each  moving  at  a 
given  uniform  velocity.  However,  the  critical  velocities  sustained 
from  this  form  of  analysis  appears  to  be  only  of  the  order  of 
magnitude  of  those  reported  in  the  literature  for  critical  velocities 
at  the  onset  of  droplet  formation.  In  reality,  the  oil  slick  layer 
has  a  finite  thickness  in  which  the  free  surface  of  the  oil  slick  does 
not  necessarily  adhere  to  the  previously  used  boundary  condition  of 
infinte  depth.  Jones  (1972)  addressed  this  problem  by  conducting  a 
linear  K-H  instability  analysis  for  a  slick  of  arbitary  but  finite, 
thickness.  Jones  concluded  that  a  thin  slick  would  be  unstable  at  a 
velocity  lower  than  the  critical  velocity  for  a  thick  slick. 

Leibovich  (1976) 


Leibovich  has  conducted  a  linear  K-H  instability  analysis 
which  leads  to  the  opposite  conclusion  from  Jones'  prediction  that  a 
thick  slick  is  more  stable  than  a  thin  one.  Leibovich  states  that 
Jones  omitted  the  effects  of  the  air-oil  surface  tension  and  air 


' 
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density,  thus  imposing  a  constant  pressure  boundary  condition  at  the 
air-oil  interface.  Leibovich  argues  that  the  justification  was  not 
valid  and  he  indicates  that  Jones  conclusions  regarding  the 
destabilizing  effects  of  decreasing  the  slick  thickness  was,  also, 
invalid.  A  brief  description  of  Leibovich's  analysis  follows. 

Considering  the  configuration  of  Figure  2.4,  the  stability 
analysis  consists  of  three  fluid  layers.  The  three  velocities  can  be 

perturbed  by  0^,  0^  and  0a*  which  must  satisfy  the  Laplace 
Equations 


a2V  =  0  a2<V  =  0  A2tfa‘  =  0 


(2.4.30) 


The  air  phase  and  water  phase  are  still  assumed  to  be 


infinite,  so  that  the  dynamic  boundary  conditions  of; 


lim  y  »  infinity  0a*  =  o  and 

lim  y  »  -infinity  0W*  _  o 


This  leads  to  the  following  solution  to  Equations(2.4.30) ; 


V  =  Aa  exp[-Ky  +  i(Kx  -  at)] 

V  =  Aw  exp[  Ky  +  i(Kx  -  at)] 


(2.4.31. a) 
(2.4.31. b) 


However,  the  finite  thickness  of  oil  depth  h  requires  a  more 


general  solution  to  the  Laplace  Equation; 


, 
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rQ  =  (A  cosh  Ky  +  Bq  sinh  Ky)  [exp  i(Kx  -at)]  (2.4.32) 

There  exists  two  interfaces  which  serve  as  streamlines 
between  the  three  liquids.  Following  a  similiar  argument  that  led  to 
the  kinematic  boundary  condition  of  Equation  (2.4.8),  the  total 
derivative  of  a  particle  path  can  be  applied  to  the  oil-air  interface 
y  -  ^2  anc*  to  the  oil -water  interface  y  =  ni ; 

Dt^  =  0  and  Dt^  =  0 
WT  JTT 

Expanding  these  derivatives  in  terms  of  the  partial 
derivatives  and  applying  them  to  both  sides  of  the  streamline 
interface  results  in  the  following  four  kinematic  boundary  conditions; 

dm'  +  U,  dm'  =  -d$'  at  y  =  h  (2.4.33. a) 

I  a  1  a 

dt  dx  dy 


dnj  + 

U0  dn{  = 

-d^ 

at  y  =  h 

(2.4. 33. b) 

dx 

dy 

dn£  + 

cz 

o 

Cl 

S3 

ro¬ 

ll 

“d^o 

at  y  =  0 

(2.4.33.C) 

dt 

dx 

dy 

dnA  + 

Uw  d„*  = 

-d  <b 1 

Ww 

at  y  =  0 

(2.4. 33- d ) 

cTt 

d7 

dy 

Two  additional  dynamic  boundary  conditions  arise  from  the 
pressure  force  balance  at  the  interfaces.  Defining  the  oil-air 
surface  tension  as  TQ^a  and  the  oil-water  surface  tension  as  TQyw, 
the  two  force  balance  equations  are; 


"  ■  '  1  f'  m 
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at  y  =  nx 


at  y  =  n2 


(2.4.34. a) 


(2.4.34.b) 


Since  the  interfaces  are  streamlines,  the  pressure  forces  can 
be  evaluated  by  a  linearized  Bernoulli  Equation,  thus,  reducing  the 
force  balance  Equation  (2.4.34. a  and  .b); 


P0(d^  +Uod*i  +9ni)-Pa(d^a  +Uad*a  +g,’l) 
dt  dx  dt  dx 

’wK,  +W  +3”2)-°o(^  +W  +gn2> 

dt  dx  dt  dx 


(2. 4. 35. a) 


(2.4.35.b) 


The  interfacial  displacements  can  be  defined  in  a  general 

form  as; 


=  a^  exp[i(Kx  -  at)]  (2.4.36. a) 

n2  =  a2  exp[i(Kx  -  at)]  (2.4.36.b) 

Substituting  the  defined  interfacial  displacements.  Equations 
(2.4.36. a  and  .b),  the  general  solutions  of  the  perturbations 
Equations  (2. 4. 31. a  and  .b)  and  (2.4.32)  into  the  six  boundary 
conditions  (Equations  (2.4.33. a  ,  .b  ,  .c  ,  and  .d)  and(2.4.35.a  and 
.b)  and  performing  the  required  differentiation  yields; 


i  (UaK  -  a)  =  AaK  exp(-Kh) 
i  al  (U0K  -  a)  =  -K(A0sinh  Kh  +  B0cosh  Kh) 


(2.4.37. a) 
(2.4.37.b) 
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i  a2  (u0K  -  a)  =  -KB0  (2.4.37.c) 

1  a2  (UwK  -  a)  =  -/^Kq  (2.4.37-d) 

poC (A0sinh  Kh  +  B0cosh  Kh)(U0K  -a)i  +ga^]  (2.4.37.e) 

-  Pa[Aai (UaK  -c)  +  gnl]  =  T0/a  aK2 

pw[Awi(UwK  -o)  +  ga2]  (2.4.37. f) 

-  poCA0i(U0K  -a)  +  ga2]  =  Ta/w  aK2 

These  six  equatons  with  the  six  unknowns  A,  A. 

B0  ,  ax  ,  and  a2  have  a  solution  only  if  the  determinant  of  the 
coefficients  vanishes.  This  leads  to  the  eigen  value  equation  for  the 
wave  speed  C  as; 

po<uo  “C )2[ Pq ( uo  -c)2tanh  Kh  +  S2(C)]  (2.4.38. a) 

+  si(C)[p0(U0  -C)2  +  S2(C)  tanh  Kh]  =  0 


Where  Sl(C)  =  0w(uw  _C)2  -  (pw-Pq)  g/K  -KTo/w 
and  S2(C)  =  0a(ua  -C)2  -  (p0-pa)  g/K  -KTo/a 

K-H  instability  occurrs  when  the  solution  for  the  wave  speed 
C  becomes  imaginary.  By  examing  the  limits  of  the  oil  thickness  h  , 
as  it  tends  to  zero  and  as  h  tends  to  infinity,  a  comparison  between 
the  relative  stability  of  thin  and  thick  oil  slicks  can  be  drawn. 

In  the  limits  of  h  tends  to  zero  (h  »  0),  Equation  (2.4.38) 


reduces  to; 
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S1(C)  +  S2(C)  =  0  or 

p  (U  -C)2  +  P  (U  -C)2  -  (p  -P  )  g  -  K(T  ,+T  ,  )  =  0 
wv  w  '  Mav  a  '  v  w  Ha;  v  o/a  o/w' 


K 

in  which  the  solution  of  the  wave  speed  C  is  (2.4.39); 


c  ■  tiV'A  *  [  1  <9  (Pw-Pa)+K(To/a+To/w))  -Pa^w,^,2]  1/2 


Pa  +  p 


w 


pa+pw  K 


p  +p 
pa  Mw 


Equation  (2.4.39)  is  the  same  solution  for  air  flow  over 
water  but  the  surface  tension  Ta/w  is  replaced  by  To/a+To/w  . 

The  form  of  the  equation  is  similiar  to  that  derived  for  an  infinite 
thickness  of  oil  (Equation  (2.4.24)).  In  the  limit  of  the  oil 
thickness  h  tending  to  zero,  the  oil  velocity  and  density  play  no  part 
in  the  stability  of  waves.  Instability  will  occur  when; 


(“wV  >  £9  +  K<To/a+To/w>]  (2’4-40> 

pa  pw  K 

The  stabilizing  influence  of  a  thin  film  on  water  as  derived 
from  the  elastic  properties  of  very  thin  films  (Lamb  Art  351)  has  been 

neglected  in  the  K-H  analysis. 

* 

Examination  of  the  limit  of  thick  oil  slicks,  when  h  tends  to 
infinity,  reduces  Equation  (2.4.38)  to; 

P0[(U0-C)2  +  S1(C)]  [p0(u0-c)2  +  S2(C)3  =  0 


(2.4.41) 


' 


. 
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When  the  first  bracket  of  Equation  (2.4.41)  is  equal  to  zero, 
the  K-H  instability  corresponding  to  the  oil-water  interface  occurs  if; 

<Uw  -  V2  >  >0  *  PW  [9  (pw  -  p0)  +  K  T0/w]  (2.4.42) 

po  pw  K 

The  minimum  critical  speed  Ucp  =  [uw  -  U0]  at  which 
instability  will  occur  is  given  for  a  minimum  value  for  the  critical 
wave  number  Kcr. 

U  2  =  2(p  +  p  )  [  T  .  g  (p  -  p  )]1/2  (2.4.43) 

cr  v  o  w  o/w  w  o 

p  p 
0  Hw 

with  a  corresponding  critical  wave  number  of; 

C  "  (»w  -  po)  9  <2-4-44> 

To/w 

The  critical  speed  for  the  thin  oil  slick  can  be  expressed  as; 

Ucr2  '  2tTo/a  *  To/w  9  fw  ]  W  (2-4‘45) 

pa  pa 

A  comparison  between  Equation  (2.4.43)  and  (2.4.45) 
indicates,  that  in  the  limits,  a  thin  oil  slick  is  considerably  more 
stable  than  a  thick  oil  slick  under  similiar  conditions.  This  result 
is  contingent  upon  the  fact  that  the  oil-air  surface  tension  has  a 
stabilizing  influence  when  the  layer  is  thin,  but  has  less  influence 
when  the  layer  is  thicker. 
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It  has  been  consistently  observed  in  various  studies  that  the 
initiation  of  oil  droplet  formation  occurs  at  velocities  in  excess  of 
those  predicted  by  K-H  analysis.  Citing  the  work  of  Drazin  (1970),  in 
which  it  was  shown  that  unstable  waves  may  be  stabilized  by  the 
non-linear  terms  ignored  in  the  linear  stability  analysis,  Liebovich 
argues  that  wave  breaking  should  occur  at  some  velocity  exceeding 

Ucr  by  the  previous  linearized  K-H  instability  analysis.  According 
to  Drazin's  non-linear  analysis,  inf initestimal  amplitude  waves  grow 

until  they  achieve  a  finite  amplitude.  The  wave  length  remains 
unchanged  during  amplification  and  corresponds  to  the  critical  wave 

number  Kcr  given  by  linear  stability  analysis  as  the  fastest  growing 
waves.  A  time  independent  equilibrium  amplitude  Ac  for  the  critical 
wave  number  Kcr  -js  described  as; 

Uw  "uo  =  Ucr (1+5/32  Kcr  2  Acr  2  )  (2.4.46) 

where  Kcr  Acr  =  s  =  wave  slope 

Liebovich  concludes  that  waves  along  the  oil-water  interface 
break  and  lead  to  droplet  formation  when  the  equilibrium,  most  stable, 
waves  achieve  a  maximum  wave  slope.  Thus  if  the  critical  breaking 

wave  slope  s  is  constant  then  droplet  formation  should  occur  at  a 

* 

velocity  as  a  constant  multiple  of  Ucr,  ie.  Uj  =  kUcr. 

Interfacial  waves  breaking  with  a  wave  slope  near  unity  would 
account  for  the  experimental  droplet  formation  velocities  reported  by 
Hydronautics.  Their  droplet  formation  velcoity  Ud  can  be  expressed 
as; 


' 


•  ,  s 
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ud  =  2.2  Ucr  (2.4.47) 

The  literature  has  repeatedly  reported  that  oil  droplet 
formation  is  basically  restricted  to  the  headwave  zone.  As  previously 
mentioned,  Wicks  has  suggested  that  the  headwave  dominance  stems  from 
the  instability  of  the  headwave  at  the  point  of  maximum  thickness 
disscussed  by  Benjamin.  Hydronautics  argue  the  droplet  formation  on 
the  lee  of  the  headwave  results  from  an  adverse  pressure  gradiant. 
Leibovich  postulates  a  more  simplistic  cause.  The  local  current  speed 
is  greatest  at  the  maximum  thickness  point.  This  locally  higher 
current  speed  implies  that  droplets  can  be  split  off  from  the  headwave 
by  breaking  of  the  K-H  waves  prior  to  oil  droplet  formation 

elsewhere.  A  range  of  Uw  _  u0  thus  exsists  where  droplets  form 
only  at  the  headwave.  At  higher  speeds  the  droplet  shedding  rate  will 
continue  to  be  greatest  at  the  headwave  but  shedding  would  be  expected 
everywhere  along  the  oil  slick.  Leibovich  derives  the  maximum  current 
velocity  from  the  recogntion  of  Benjamin's  headwave  Froude  Number  as 
Umax  =  2Ud.  The  previous  section  has  shown  that  the  current 
velocity  for  droplet  formation  is  Ud  =  2.2  Ucr,  thus  the  actual 
velocity  at  which  the  droplets  are  torn  from  the  oil  slick  is  given  as 

Un,ax  =  2.2  (2)1/2  Ucr  =  3.1  Ucr  (2.4.48) 

Leibovich  concludes  that  since  droplets  will  not  be  shed  from 
the  downstream  portion  of  the  oil  slick  until  the  local  current 
velocity  exceeds  Umax>  there  is  a  range  of  currents  velocities  Ud; 

2*2  Uc^  <  Ud  <  3.1  Ucr 


' 
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for  which  droplet  formation  must  be  confined  to  the  headwave. 


Milgram  and  Van  Houten  (1978) 


Milgram  and  Van  Houten  have  reported  several  experimental 
observations  that  reinforce  the  hypothesis  that  droplets  are  generated 
by  the  breaking  of  Kelvin  -  Holmholtz  waves.  Initially,  they  were 
concerned  that  the  observed  interfacial  waves  were  results  of  the 
Tol Imien-Schl ichting  instability  phenomena.  Tol lmien-Schl ichting 
instability  occurs  during  the  transition  from  laminar  to  turbulent 
boundary  layer  flow  along  the  interface  when  the  viscous  forces,  which 
have  been  restraining  the  small  disturbances,  are  overwhelmed  by  the 
intertial  dynamic  forces.  Milgram  and  Van  Houten  present  the 
following  three  reasons  to  indicate  that  the  observed  interfacial 
waves  were  in  fact,  Kel vin-Helmhol tz  waves  and  not  the  results  of 
fluctuations  in  Tollmien-Schl ichting  waves. 


1)  The  characteristics  of  Kelvin-Helmholtz  waves  depend  on 
fluid  densities,  interfacial  tensions,  and  flow  speed,  but 
not  on  downstream  location.  For  a  given  oil  type  and  flow 
speed,  the  waves  we  observed  had  constant  phase  speed  and 
wavelength  as  they  moved  downstream,  thus  behaving  as 
Kelvin-Helmholtz  waves.  On  the  other  hand,  Tollmien- 
Schl  ichting  waves  would  be  expected  to  become  longer  and 
move  more  slowly  as  they  move  downstream  as  a  result  of  the 
increasing  Reynolds  Number  and  boundary '1 ayer  thickness  with 
downstream  position. 

2)  As  the  flow  speed  was  slowly  increased,  the  waves  first 
appeared  near  the  leading  edge  of  the  oil  slick,  as  would  be 
expected  for  Kelvin-Helmholtz  waves  because  the  shear  layer 
in  the  water  is  stronger  there  than  it  is  further 
downstream.  Tol Imein-Schl ichting  waves  would  be  expected  to 
appear  downstream  first  owing  to  the  larger  Reynolds  Number 
based  on  slick  length  for  downstream  regions  than  for 
upstream  regions. 


I 
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3)  At  the  higher  flow  speeds  when  the  interface  is  quite 
rough,  starting  a  short  distance  downstream  from  the  region 
where  the  waves  begin,  the  boundary  layer  must  become 
turbulent  and  cease  to  support  Tollmien-Schl ichting  waves 
having  any  regularity.  However,  we  observed  interfacial 
waves  possessing  some  regularities  over  almost  the  entire 
slick,  although  their  amplitudes  were  substantially 
diminished  far  downstream.  This  indicates  that  they  were  not 
Tol lmien-Schl ichting  waves. 


Milgram  and  Van  Houten  observed  that  the  waves  did  not  appear 
at  the  leading  edge,  but  rather  appeared  a  short  distance  aft  of  the 
leading  edge.  They  attribute  this  to  the  stabilizing  effect  of  a  thin 
slick  as  indicated  by  Leibovich,  since  the  leading  edge  of  the  oil 
slick  is  very  thin  in  comparision  to  the  headwave  maximum  thickness. 
The  amplitude  of  the  interfacial  waves  was  largest  in  the  vicinity  of 
the  headwave  corresponding  to  the  region  of  maximum  shear  coefficients 
(section  2.2.)  and  both  shear  and  wave  amplidute  was  found  to  increase 
with  increasing  water  velocity. 


Droplet  formation  always  occured  on  the  leeward  side  of  the 
headwave  where  the  local  water  velocity  was  highest.  However,  under 
certain  conditions  some  of  the  oil  slicks  had  a  geometry  with  the 
highest  local  velocity  occurring  in  the  viscous  zone  of  the  oil  slick 
although  K-H  waves  did  not  break  within  the  viscous  zone.  From  this 
observation,  Milgram  and  Van  Houten  suggest  that  the  droplet 
generation  phenomena  is  more  complex  than  postulated  by  Leibovich,  who 
determined  that  entrainment  ocurrs  with  any  specific  oil  wherever  a 
certain  critical  speed  is  exceeded  in  the  flow  beneath  the  oil, 
independent  of  the  portion  of  the  slick  at  which  the  required  velocity 
occurs.  Milgram  and  Van  Houten  speculate  that  the  amplitude  of  the 


■ 
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interfacial  waves  depends  not  only  on  the  water  speed  below  the 
boundary  layer  in  the  water,  but  also  on  the  boundary  layer  thickness 
itself.  When  the  layer  is  thick,  changes  in  interfacial  shape  lead 
to  smaller  pressure  changes  than  when  the  boundary  layer  is  thin. 
Droplet  formation  was  observed  only  in  the  headwave  where  the  local 
water  velocity  was  high  and  the  boundary  layer  was  thin  (as  indicated 
by  high  friction  coefficients).  Downstream,  though  the  local  velocity 
was  high,  the  boundary  layer  was  considerably  thicker  and  no  droplet 
formation  was  observed.  The  wave  slope  of  breaking  interfacial  waves 
was  reported  to  be  close  to  unity  as  required  by  Leibovich's  theory, 
but  increased  boundary  layer  thickness  reduces  the  amplitude  growth 
more  than  was  suggested  by  Leibovich. 

Another  aspect  of  K-H  instability  reported  by  Mi  gram  and  Van 
Houten  was  the  apparent  inconsistency  between  the  most  unstable  wave 
length  of  classical  K-H  waves  and  those  observed  during  the  course  of 
their  experiments.  They  cite  observing  wave  lengths  in  the  range  of 
10  to  25  centimeters  for  a  heavy  mineral  oil,  whereas,  the  critical 
wave  length  given  by  Equation  (2.4.43)  indicates  the  waves  should  have 
been  4.3  centimeters  in  length.  Evidently  the  relatively  large 
damping  of  the  very  short  waves  prevented  them  from  growing  large 

* 

enough  to  break.  Considerable  shorter  waves  were  observed  in  diesel 
oil,  in  which  Equation  (2.4.43)  yielded  a  critical  wave  length  of  2.4 
centimeters.  Since  the  kinematic  viscosity  of  the  diesel  oil  is  less 
than  3  percent  of  that  for  the  heavy  mineral  oil,  the  shorter  K-H 
waves  undergo  less  damping  and  as  a  result  waves  of  the  order  of  Kcr 
were  observed.  Milgram  and  Van  Houten  concluded  that  the  viscosity  of 


■ 
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the  oil  plays  a  role  in  determing  the  critical  speed  of  droplet 
formation  by  suppressing  the  growth  of  the  most  unstable  K-H  waves 
when  the  viscosity  is  sufficiently  high. 

2.5  Summary 


In  summary,  the  following  conclusions  pertaining  to  the 
headwave  zone  can  be  drawn.  The  theoretical  analysis  by  Benjamin 
appears  to  describe  the  asymptomatic  thickness  of  the  neck,  downstream 
of  the  maximum  thickness  of  the  headwave,  with  a  fair  degree  of 
accuracy.  This  analysis,  based  on  the  concept  of  momentum  balance, 
allows  for  the  prediction  of  the  headwave  with  a  simple  Froude  number 
model.  Benjamin's  analysis  further  describes  the  magnitude  of  the 
energy  loss  which  is  required  for  this  momentum  balance.  However,  the 
neglect  of  shear  stresses  results  in  the  breaking  of  the  headwave  to 
achieve  this  desired  energy  dissipation.  In  order  for  the  headwave  to 
break,  the  oil  thickens  to  roughly  twice  the  downstream  depth  causing 
the  water  to  become  super-critical  locally.  The  results  of  Milgram 
and  Van  Houten's  testing  strongly  suggests  that  the  shear  stresses 
predominate  in  the  energy  loss  mechanism  and,  thus,  at  the  lower 
velocities,  and  consequently  lower  energy  loss  levels,  the  oil  slick 
is  observed  without  a  breaking  headwave.  However,  at  the  larger 
velocities  when  the  energy  loss  requirements  are  maximum,  the  shear 
stresses  are  insufficient  to  consume  all  the  energy  and  the  headwave 
must  break  in  the  manner  previously  described. 
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There  is  total  agreement,  between  the  various  investigators, 
concerning  the  governing  forces  which  come  into  play  within  the 
viscous  zone.  There  are  several  models  developed  which  describe  the 
growth  of  the  viscous  zone  by  various  methods  of  accumulating  the 
interfacial  shear  stresses.  It  is  only  magnitude  of  shear  stress  and 
the  method  of  utilizing  the  shear  stresses  that  differs. 

The  simplest  models  employing  constant  shear  were  proposed  by 
Cross  and  Hoult,  by  Lindenmuth,  Miller,  and  Hsu,  and  by  Hale,  Norton, 
and  Rodenberger.  This  last  study  included  serveral  empirical 
relationships  to  describe  the  various  aspects  of  the  viscous  zone. 

There  are  two  main  attributes  which  stem  from  Wicks' 
analysis.  The  evaluation  of  interfacial  shear  stress  is  not  limited 
to  a  spacial  constant  coefficient  of  friction.  Furthermore,  the  shear 
stress  is  evaluated  from  the  difference  between  the  interfacial 

velocity  and  the  free  stream  velocity  U.  The  major  contribution 
of  Wicks'  analysis  has  been  the  attempt  to  include  the  mechanics  of 
flow  within  the  oil  slick  in  evaluation  of  the  oil  slick  thickness  in 
the  viscous  zone. 

The  principal  contribution  of  Wilkinson's  viscous  zone  model 
is  the  employment  of  the  boundary  shear  for  the  flow  of  finite  depth. 
The  major  practical  application  of  this  concept  should  be  in  the 
extrapolation  of  model  to  prototype  oil  slick  geometries. 
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Milgram  and  Van  Houten's  detailed  analysis  of  the  viscous 
zone  demonstrated  conclusively  the  importance  of  the  interfacial  shear 
stress  in  the  development  of  the  viscous  zone.  They  suggested  that 
the  friction  distributions  fall  into  two  distinct  types  but  were 
unable  to  formulate  boundary  limits. 

The  onset  of  entrainment  failure  is  the  result  of  Kelvin  - 
Helmholtz  instability  at  the  oil-water  interface.  The  most 
comprehensive  analysis  of  linear  K-H  instability  is  given  by 
Leibovich.  Linear  K-H  instability  anaylsis  under  predicts  the 
critical  velocity  for  interfacial  breakup.  Citing  the  work  of  Drazin 
(non-linear  K-H  instabilities),  Leibovich  demonstrates  the  interfacial 
breakup  occurs  at  a  finite  value  higher  than  predicted  by  linear  K-H 
analysis.  Milgram  and  Van  Houten  further  argue  that  the  development 
of  a  interfacial  boundary  layer  retards  the  critical  breakup  velocity 
by  suppressing  the  growth  of  the  shortest  (most  unstable)  waves. 

This  literature  review  documents  the  basic  oil  slick 
phenomenon  for  the  ice  free  oil  slick.  The  subsequent  chapters  will 
examine  various  components  of  the  oil  slick  as  they  are  affected  by 
the  addition  of  a  solid  ice  cover. 
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CHAPTER  3 


SCOPE  OF  LABORATORY  STUDY 


3.1  Introduction 


The  experimental  phase  of  this  study  was  directed  to  the 
investigation  of  the  characteristic  behaviour  of  ice  covered  oil 
slicks.  The  experiments  were  designed  to  address  the  following 
objectives: 

1)  Initially,  to  determine  insight  into  the  basic 
characteristics  of  oil  slicks,  control  experiments  were 
conducted  under  a  varying  range  of  flow  conditions  in  order 
to  observe  and  document  the  behaviour  of  ice  free  oil  slicks. 

2)  A  primary  objective  was  to  identify  the  fundamental 
characteristics  of  oil  slicks  as  influenced  by  the  presence 
of  an  ice  cover.  Under  a  similar  set  of  flow  conditions 
utilized  in  the  control  experiments,  the  addition  of  a  model 
ice  cover  permitted  the  observation  and  documentation  of  the 
basic  differences  in  behaviour  of  ice  covered  oil  slicks. 

3)  In  order  to  determine  the  mechanisms  by  which  the 
differences  in  behaviour  evolved,  detailed  velocity  profiles 
within  the  oil  slick  and  water  phase  were  measured  in  both 
the  control  and  ice  covered  experiments.  This  provided  the 
basis  for  an  evaluation  of  the  shear  stress  at  the  interface 
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and  examination  of  its  changes  as  influenced  by  an  ice  cover. 

4)  The  formation  of  interfacial  waves  has  been  demonstrated 
to  be  responsible  for  the  initiation  of  entrainment  failure 
of  oil  slicks.  An  objective  in  this  study  was  to  indentify 
the  role  of  a  solid  boundary  (ice  cover)  as  it  affects  growth 
of  interfacial  waves. 

In  order  to  pursue  these  objectives,  several  measurements 
were  conducted  to  establish  the  various  properties  of  oil  slick 
behaviour  relating  to  the  slick  thickness  profile,  velocity  profile, 
and  the  interfacial  wave  characteristics.  Both  the  control  and  ice 
covered  experiments  encompassed  a  complete  range  of  Froude  Numbers  in 
which  the  primary  variable  was  the  water  current. 

The  remainder  of  the  chapter  describes  the  various  aspect  of 
the  testing  arrangement  and  outlines  the  experimental  procedures 
utilized  in  this  study. 

3.2  Equipment 


The  experimental  tests  were  conducted  at  the  T  Blench 
Hydraulics  Laboratory,  located  at  the  University  of  Alberta. 

The  model  basin  utilized  for  all  the  oil  slick  tests 
consisted  of  a  self  contained  re-circulating  flume.  The  useable  test 
section  measures  5.0  m  in  length,  0.46  m  inside  width,  and  0.9  m  in 
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depth.  Both  of  the  flume  walls  were  constructed  of  plexiglass,  12  mm 
in  thickness,  which  facilitated  both  ease  of  observation  and 
photography.  The  bed  of  the  flume  consisted  of  aluminum  plate 
constructed  in  a  fixed  postion  of  zero  slope.  A  schematic  diagram  of 
the  general  flume  configuration  is  depicted  in  Figure  3.1. 

Located  at  the  upstream  end  of  the  test  section  was  a  headbox 
measuring  approximately  1.8  m  in  length,  1.0  m  in  depth  and  0.5  m  in 
width.  Water  was  introduced  from  a  pipe,  over  the  full  width  of  the 
bottom,  at  the  upstream  portion  of  the  headbox.  Inside  the  headbox 
was  a  baffle  system  consisting  of  a  full  depth  coarse  metal  screen 
covered  with  a  rubberized  porous  material  similar  to  a  furnace  air 
filter.  This  baffle  system  coupled  with  a  floating  plywood  board  was 
required  to  minimize  the  severe  turbulent  fluctuation  of  the  water 
surface  induced  by  the  water  discharging  into  the  headbox. 

The  sufficiently  narrow  width  of  the  headbox  removed  the  need 
for  a  laterial  transition.  Between  the  headbox  and  the  upstream  of 
the  test  section,  the  flume  was  fitted  with  an  entrance  gate.  During 
the  course  of  testing,  this  gate  was  raised  above  the  water  surface 
and  consequently  had  no  effect  on  the  control  of  the  flow.  Water 
levels  in  the  flume  were  maintained  by  the  adjustment  of  the  variable 
position  undershot  gate,  located  downstream  of  the  test  section. 

Water  discharging  from  the  undershot  gate  was  conveyed  down  a  short 
free  flowing  chute  into  the  sump. 

The  steel  sump,  isolated  solely  to  the  flume,  provided  a 
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storage  capacity  of  4.5  m^.  The  self  contained  re-circulating  flume 
was  furnished  with  a  pump  and  appropriate  piping  between  the  sump  and 
headbox.  The  pump,  in  its  current  piping  configuration  was  capable  of 
delivering  a  maximum  discharge  of  0.061  m3/s.  The  flow  was 
regulated  by  a  gate  valve  downstream  of  the  pump  adjacent  to  the 
headbox.  The  regulated  discharge  was  monitored  with  an  inline 
magnetic  flowmeter  located  between  the  pump  and  the  gate  valve. 

The  afore-mentioned  flume  configuration  had  been  constructed 
several  years  prior  to  the  initiation  of  this  study.  This  flume  was 
selected  for  use  in  this  study,  primarily  because  its  isolated  sump 
removed  the  danger  of  lost  oil  contaminating  the  rest  of  the 
hydraulics  laboratory.  However,  in  order  to  conduct  the  oil  spill 
tests,  two  additional  features  were  incorporated  into  the  flume 
configuration. 

The  first  modification  consisted  of  constructing  an  auxiliary 
undershot  gate  within  the  test  section.  Located  1.0  m  upstream  of  the 
water  level  control  gate,  this  second  gate  provided  an  oil  slick  boom 
of  variable  boom  draft.  A  plywood  sheet  was  mounted  on  a  rack  and 
pinion  thus  allowing  adjustment  in  the  boom  draft.  Since  the  boom  was 
free  to  move  in  the  vertical  direction,  rubber  strips  3  cm  wide,  were 
mounted  at  the  side  of  the  boom  on  the  upstream  face  forming  a  concave 
seal  against  the  wall  of  the  flume. 

The  water  flow  in  a  flume  has  an  inherent  lateral  velocity 
distribution  resulting  from  the  boundary  layer  effects  of  the  flume 
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walls.  This  velocity  pattern  yields  a  stagnation  pressure 
distribution  at  the  boom  associated  with  a  lateral  vortex  formation  at 
the  walls  of  the  flume  immediately  upstream  of  the  boom.  If 
unarrested,  these  vortices  quickly  drain  the  oil  slick  beneath  the 
boom,  regardless  of  the  boom  draft.  Suppression  of  these  vortices  was 
accomplished  by  the  installation  of  a  0.1  m  thick  strip  of  rubberized 
porous  mat  on  the  upstream  face  of  the  boom.  This  mat  reduced  the 
higher  water  velocities  over  the  central  portion  of  the  boom  and 
suppressed  the  resulting  differential  pressure  gradient  to  a  level 
incapable  of  driving  the  vorticies.  A  secondary  functon  of  the  porous 
mat  was  to  absorb  wave  energy  and  minimize  the  reflection  of  waves 
(characteri Stic  of  a  rigid  boom)  in  the  oil. 

The  second  addition  to  the  flume  was  the  incorporation  of  a 
removable  ice  cover.  The  model  ice  cover  consisted  of  a  single  0.2  m 
x  0.012  m  thick  strip  of  plexiglass  and  two  strips  of  plywood  0.13  m 
x0.012  m  thick  fastened  on  either  side  of  the  plexiglass  to  form  a 
solid  smooth  cover  0.46  m  wide  x  0.012  m  thick.  The  central 
transparent  section  of  the  ice  cover  permitted  both  the  observation 
and  the  overhead  lighting  of  the  oil  slick.  The  model  cover  was  3.5  m 
long  and  was  floated  on  the  surface  of  the  oil  and  water  upstream  of 
the  boom. 


The  model  cover  was  narrower  than  the  flume  to  facilitate 
ease  of  removal,  thus  the  ice  cover  was  propelled  upstream  by  the 
shear  forces  exerted  by  the  oil  slick  as  the  upper  layer  of  oil  flowed 
upstream.  The  ice  cover  was  wedged  in  place  to  prevent  this  problem. 
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Since  the  model  ice  cover  consisted  of  only  a  smooth  solid 
sheet  it  did  not  attempt  to  model  all  aspects  of  a  real  ice  cover. 

Only  those  characteristics  determined  by  a  solid  boundary  on  the 
surface  of  the  oil  were  investigated  and  no  consideration  was  given  to 
behavioural  changes  induced  by  the  interaction  of  the  oil  and  ice  with 
respect  to  their  thermal  and  or  chemical  properties. 

3.3  Measurements 


To  achieve  the  experimental  objectives  outlined  at  the 
begining  of  this  chapter,  three  distinict  characteristics  of  the 
oil-slick  were  monitored.  The  principal  measurement  consisted  of 
determining  the  chartacteristic  length  and  thickness  profile  over  a 
complete  range  of  flow  conditions. 

A  major  emphasis  of  this  study  addressed  the  measurement  of 
velocity  profiles  both  in  and  below  the  oil  slick.  This  type  of 
measurement  was  fundamental  in  establishing  the  interfacial  shear 
stresses  and  in  developing  an  understanding  of  the  internal  flow 
mechanisms  of  an  oil  slick. 

t 

The  last  series  of  measurements  consisted  of  defining  the 
interfacial  wave  characteristics  as  functions  of  both  time  and  space. 

The  remainder  of  the  section  addresses  the  specific  details 
involved  in  the  three  forms  of  measurements. 
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The  characteristic  length  scales  of  the  oil  slick  were 
successfully  determined  by  two  independant  methods.  The  most  obvious 
of  these  was  to  measure  the  oil  slick  profile  outlined  on  the 
plexiglass  flume  walls.  The  oil  slick  was  in  continous  motion  as 
interfacial  waves  travelled  down  its  entire  length.  Thus,  the 
thickness  profile  was  recorded  by  visual  integration  (with  time)  of 
the  interfacial  position  every  0.05  m  along  its  length.  Although  the 
oil  slick  shape  at  the  flume  walls  appeared  to  be  adequately 
representative  of  the  oil  thickness  for  the  full  width  at  a  given 
section,  the  leading  edge  of  the  slick  was  not  spacial  invariant.  The 
recorded  slick  length  was  derived  from  the  visual  averaging  the 
positon  of  the  leading  edge  of  the  oil  slick  over  the  flume  width. 

As  second  set  of  oil  slick  profiles  was  determined, 
indirectly  by  recording  the  oil  slick  phenomena  on  super  8  mm  film. 

To  facilitate  accurate  measurements,  the  flume  wall  was  outlined  with 
thin  black  tape  to  form  a  grid  network  composed  of  0.05  m  squares. 
Viewing  of  the  film  enabled  a  frame  by  frame  analysis  of  the  oil  slick 
profiles,  and  facilitated  describing  the  instantaneous  interfaciial 
wave  forms. 

* 

A  third  method  employed  to  measure  the  oil  slick  profile  was 
somewhat  unproductive.  This  procedure  consisted  of  utilizing  a  water 
surface  probe  to  indicate  the  oil-water  interface.  The  probe 
consisted  of  a  thin  steel  rod  in  which  the  resistance  of  water  at 
varying  depths  of  penetration  provided  input  for  a  feed-back  circuit 
which  in  turn  positioned  the  probe  mechanically  to  maintain  a  constant 
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depth  of  penetration.  Since  oil  is  a  relative  insulator,  the 
electronic  controls  could  be  adjusted  to  have  the  probe  penetrate  the 
oil  and  detect  the  oil-water  interface.  This  type  of  apparatus  could 
ideally  produce  a  continuous  plot  of  the  oil-water  interface  in  the 
middle  of  the  flume  away  from  any  wall  effects.  However,  with  the 
passage  of  an  interfacial  wave,  the  probe  would  follow  the  water 
surface  down  as  the  oil  thickened  but  it  failed  to  return  as  the  wave 
passed.  The  apparent  reason  for  the  probe's  failure  was  that  a  thin 
layer  of  oil  would  adhere  to  the  probe  and  fail  to  rise  at  the  same 
rate  as  the  general  oil-water  surface  owing  to  the  oil's  high 
viscosity.  Thus  with  a  coating  of  oil  on  the  probe,  the  probe 
constantly  stayed  at  the  maximum  oil  depth  with  successive  passage  of 
interfacial  waves 

One  of  the  principal  objectives  of  this  study  was  to  examine 
the  interfacial  shear  stresses  as  influenced  by  the  addition  of  an  ice 
cover  to  the  oil  slick.  Earlier  investigators  have  evaluated  the 
interfacial  shear  stress  indirectly  from  an  overall  force  balance  of 
the  oil  slick.  Wicks  (1968)  considered  the  effects  of  a  velocity 
distribution  by  assuming  a  linear  distribution  throughout  the  oil 
slick.  Wicks'  velocity  profile  consisted  of  identical  velocities  and 
quantities  of  flow  above  and  below  the  midpoint  of  the  oil  slick. 
Further  researchers  have  shown  the  velocity  distribution  to  deviate 
drastically  from  Wicks'  simplified  model,  but  their  velocity 
measurements  have  consisted  only  of  crude  estimates  of  the  top  surface 
and  interfacial  surface  velocities.  Therefore  the  major  emphasis  of 
this  study  is  to  delinate  the  internal  velocity  distribution  of  an  oil 
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slick  and  evaluate  its  role  on  the  interfacial  shear  stress. 

In  this  investigation,  three  independent  techniques  were  employed 
to  describe  the  velocity  field  within  the  oil  slick.  Briefly  they 
were: 


1)  The  timing  of  particles  on  both  the  top  surface  and 
interfacial  surface  was  utilizied  to  determine  each  of  the 
velocities  respectfully. 

2)  A  pitot  tube  coupled  with  a  differential  pressure 
transducer  was  used  to  establish  the  velocity  field  below  the 
oil -water  interface. 

3)  A  laser  anemometer  was  used  to  accurately  establish 
velocity  profiles  in  the  water  phase  and  within  the  oil  slick 
itself. 

These  first  two  methods  are  similar  to  those  utilized  by  the 
previous  researchers.  However,  due  to  the  major  drawbacks  inherent  in 
their  useage  they  are  limited  to  the  crude  estimation  of  the  velocity 
distribution.  The  timing  of  particles  has  the  obvious  limitation  of 
establishing  only  two  points  on  a  velocity  distribution.  The  major 
handicap  associated  with  the  use  of  a  pitot  tube  is  that  the  size  of 
the  instrument  is  so  large,  compared  to  the  flow,  that  the  instrument 
invariably  interacts  with  and  induces  changes  in  the  velocity  field. 

On  the  other  hand,  a  laser  anemonter  does  not  disturb  the 
flow  and  can  be  used  to  accurately  define  the  velocity  of  a  point.  It 
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was  preceived  that  the  use  of  a  laser  anemonter  would  be  unique  to 
this  oil  slick  study  and  it  would  provide  considerable  data  on  the 
velocity  profile  in  the  oil  slick. 

A  brief  description  of  each  of  these  three  methods  follows. 

The  first  method  of  measuring  the  travel  time  of  particles 
was  utilized  in  most  oil  slick  tests,  including  those  not  directed 
solely  to  velocity  measurements.  This  method,  though  somewhat  crude, 
provided  a  quick  and  easy  verification  of  the  velocities  determined  by 
the  other  two  methods. 

Each  velocity  measurement  at  a  given  point  consisted  of  the 
average  of  a  minimum  of  three  separate  trials.  An  individual  trial 
involved  the  release  of  a  small  paper  flag  upstream  to  the  point  of 
interest  and  timing  its  travel  over  a  distance  of  at  least  30  cm.  The 
corresponding  velocity  was  assigned  to  the  midpoint  for  the  test  reach. 

The  small  paper  flag  would  float  on  the  water  at  the  oil- 
water  interface.  However,  in  some  cases,  the  paper  flags  released  at 

the  interface  would  assume  an  orientation  deviating  from  parallel  with 

* 

the  interface.  This  behaviour  would  result  in  an  observed  velocity 
slightly  lower  than  the  true  interfacial  velocity.  This  problem  was 
not  evident  with  the  top  surface  and  it  was  considered  that  the 
surface  flags  accurately  represented  the  surface  velocity. 


This  method  was  subsequently  used  to  estimate  point 
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velocities  for  both  the  interface  and  top  surface  at  intervals 
approximately  60  cm.  apart. 

The  pitot  tube  measurements  were  orginally  used  as  an 
independent  check  of  the  laser  anemometer  velocity  measurements 
conducted  below  the  oil  in  the  water  phase.  However,  an  attempt  to 
measure  the  velocities  using  a  pitot  tube,  within  the  oil  slick  was 
conducted  with  some  degree  of  success. 

The  pitot  tube  was  connected  to  a  differential  pressure 
transducer  with  a  full  scale  range  of  ±  1.0  inches  of  water.  The 
analog  output  from  the  transducer  was  input  to  a  strip  chart  recorder. 
Individual  point  velocities  were  derived  from  averaging  the  strip 
chart  recording  of  a  two  or  three  minute  continuous  sample. 

There  were  two  major  problems  associated  with  the  utilization 
of  the  pitot  tube  measurements  in  this  study.  First,  the  extension  of 
the  pitot  tube  into  the  water  from  the  top  surface  required  that  the 
tube  penetrate  through  the  oil  slick  into  the  water  phase.  The 
portion  of  the  tube  penetrating  onto  the  water  below  the  oil 

experienced  a  form  drag  associated  with  the  water  velocity.  This  form 

* 

drag  would  result  in  a  reduction  in  pressure  within  the  immediate 
vicinity  of  the  pitot  tube.  Consequently  some  of  the  oil  from  the 
slick  would  be  "sucked"  into  the  water  phase.  However,  this  problem 
was  restricted  to  a  finite  layer  below  the  oil  slick.  Below  this 
layer  the  intrusion  of  oil  was  insignificant,  compared  to  the  depth  of 
velocity  measurement,  so  the  oil  failed  to  adversely  influnce  the 
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reading.  Above  this  finite  layer,  in  close  proximity  to  the 
interface,  the  pitot  tube  had  insufficient  penetration  resulting  in  a 
minimum  form  drag  and  no  noticeable  disturbance  in  the  oil-water 
interface. 

The  second  major  drawback  concerning  the  usage  of  pitot  tube 
occurred  during  the  measurements  of  the  interfacial  velocities.  The 
interface  exists  in  a  continously  changing  form  in  the  presence  of 
interfacial  waves,  travelling  the  entire  length  of  the  oil  slick. 

With  time  this  results  in  a  change  in  the  fluid  medium  at  a  given 
fixed  depth  of  penetration  in  the  vicinity  of  the  interface.  A 
problem  resulted  with  the  passage  of  interfacial  waves  as  the  static 
and  dynamic  ports  of  the  pitot  tube  were  exposed  to  either  oil  or 
water.  The  oil  and  water  phases  would  transmit  pressure  fluctations 
as  different  rates  resulting  in  some  erroneous  readings  at  the 
interface. 

Even  with  these  major  obstacles  inherent  with  the  usage  of  a 
pitot  tube  at  the  interface,  the  measurements  proved  to  confirm  the 
measurements  obtained  by  the  laser  anemometer. 

* 

The  laser  dopier  anemometer  (LDA)  is  a  sophisticated  and 
accurate  technique  of  measuring  fluid  velocities.  Even  though  it  is  a 
complex  and  extremely  expensive  method,  it  offers  the  primary 
advantage  of  non-disturbance  within  the  flow,  since  only  the  laser 
beam  is  in  direct  contact  with  the  flow.  The  flow  in  an  oil  slick  is 
confined  to  a  small  layer  and  any  conventional  probe  induces  changes 
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in  the  velocity  field  as  it  attempts  to  measure  it.  It  is  the 
non-disturbance  property  of  laser  anemometery  that  makes  it 
attractive  for  this  type  of  study. 

A  second  advantage  to  usage  of  this  technique  is  its  possible 
utilization  when  the  condition  of  a  solid  ice  cover  exists  on  top. 

The  equipment  used  in  laser  anemometery  consists  of  a  laser, 
transmitting  optics,  receiving  optics,  and  a  signal  processor.  A 
simple  model,  termed  the  fringe  model,  describes  the  basic  physics 
utilized  in  laser  anemometery. 

First  the  laser  beam  is  split  with  a  prism  into  two  parallel 
beams.  These  beams  are  focused  with  a  diverging  lens  characterized  by 
a  short  focal  length,  to  form  an  intersection  point  within  the 
velocity  field.  If  a  projection  of  the  intersection  were  conducted, 
the  image  formed  would  be  a  pattern  of  vertical  lines  of  light  and 
dark.  These  are  the  interference  fringes  caused  by  the  intersection 
between  the  two  monochromic  laser  beams.  The  fringe  spacing  is  a 
function  of  the  intersection  angle  and  the  wave  length  of  the  laser 

light.  As  a  particle  within  the  velocity  field  passes  through  the 

* 

intersection  point  it  reflects  scattered  light  as  it  moves  from  the 
dark  to  light  fringes.  This  "signal"  is  recieved  by  an 
electro-optical  device  and  is  converted  to  an  analog  electrical 
signal.  The  signal  takes  the  form  of  a  sine  wave  with  an  amplitude 
proportional  to  the  amount  of  light  scattered  in  a  period  equal  to  the 
time  taken  by  the  particle  to  cross  a  single  fringe  of  light  and 
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dark.  Since  the  fringe  spacing  is  constant  with  a  given  intersection 
angle*  the  instantaneous  velocity  can  be  calculated  from  the  signal 
period. 


There  are  several  methods  of  utilizing  the  laser  anemometery 
technique.  The  first,  termed  forward  scatter,  exists  when  the 
transmitting  and  receiving  optics  are  on  opposite  sides  of  the  flume. 
The  second  method,  termed  back  scatter,  exists  when  both  the 
transmitting  and  receiving  optics  reside  on  the  same  side  of  the 
flume.  The  forward  scatter  method  is  preferable  since  the  resulting 
signal  is  much  stronger  and  easier  to  track.  However,  physical  space 
limitation  on  the  flume  dictated  the  use  of  the  back  scatter  method. 

As  a  result  of  using  back  scatter,  and  the  resulting  weaker  signal, 
the  maximum  penetration  of  the  intersection  point  was  limited  to 
approximatly  5  cm  inside  of  the  flume  wall.  This  could  have  been 
easily  increased  with  the  use  of  a  stronger  laser  had  one  been 
available.  One  advantage  gained  with  the  back  scatter  method  is  that 
it  is  fairly  simple  to  mount  both  the  transmitting  and  receiving 
optics  on  the  same  transversing  base. 

In  addition  to  the  limitation  of  only  5  cm  penetration  of  the 

t 

oil  slick  for  the  velocity  measurements,  a  second  handicap  was 
encountered  when  measuring  zero  velocity  in  the  oil  slick  layer. 
Interfacial  waves  travelling  the  length  of  the  oil  slick  induce  a 
superimposition  of  small  orbital  velocities  over  the  shear  velocities 
of  the  oil  slick  in  its  steady  state.  Therefore,  in  certain  regions 
the  instantaneous  velocity  at  a  given  point  could  be  either  upstream 
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or  downstream.  In  the  laser  anemometer  configuration  utilized  during 
this  study,  the  instrumentation  was  unable  to  detect  the  direction  of 
velocity.  Rather  the  calculations  determined  were  limited  to  the 
absolute  magnitude  of  velocity.  This  was  of  minor  consequence  when 
the  oil  velocity  direction  was  obvious  but  in  regions  of  near  zero 
velocity,  when  the  instantaneous  velocity  alternated  direction,  the 
time  average  velocity  at  a  given  point  consisted  of  a  summation  of  the 
magnitude  of  instantaneous  velocities.  Thus,  under  certain  conditions 
when  the  time  average  velocity  was  zero  resulting  from  the 
cancellation  of  the  upstream  and  downstream  component  velocities,  the 
LDA  recorded  erroneous  velocity  magnitude. 

A  method  exists  to  overcome  this  handicap  but  the  required 
instrumentation  was  unavailable  for  the  study.  Briefly,  it  consists 
of  the  measurement  of  velocity  deviations  from  an  artifical  velocity 
base  greater  than  any  encountered  in  the  velocity  field. 

Though  some  problems  were  encountered  in  the  utilization  of 
the  laser  anemometer,  the  final  results  proved  to  meet  expectations. 
The  use  of  a  LDA  system  allows  for  the  definition  of  the  time  average 
velocity  gradient  encountered  in  the  internal  flow  of  oil.  This 
technique  was,  also,  implemented  when  an  ice  cover  extended  the  oil 
slick.  Chapter  4  presents  the  experimental  velocity  profiles 
depicting  the  difference  in  oil  slick  circulation  as  influenced  by  the 
addition  of  an  ice  cover. 

The  investigation  of  interfacial  waves  on  the  oil  slick  was 
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composed  of  two  forms  of  measurements.  First,  wave  speeds  were 
estimated  from  the  timing  of  wave  crests  passing  fixed  points.  A 
deficiency  in  these  measurements  stemmed  from  the  continuously 
changing  form  of  the  waves.  The  interfacial  waves  expanded  and 
diminished  in  amplitude  as  they  migrated  downstream.  Only  the  largest 
waves  would  maintain  a  recognizable  shape  for  a  continuous  distance, 
and  thus  only  the  largest  waves  could  be  monitored. 

A  second  and  more  detailed  analysis  of  the  interfacial  waves 
was  possible  with  the  use  of  super  8  mm  film.  The  side  of  the  flume 
was  marked  out  into  a  grid  pattern.  Knowledge  of  the  filming  speed 
thus  allowed  for  a  detailed  frame  by  frame  examination  of  the  film 
footage.  These  measurements  resulted  in  both  definition  of  wave 
velocity  and  characteristic  shape. 

3.4  Oil  Supply 


The  basic  characteristics  of  oil  slick  behaviour  vary 
according  to  properties  of  the  particular  oil  involved.  Three 
distinct  types  of  oil  were  utilized  in  the  experimental  phase  of  this 
study.  The  physical  properties  of  these  oil,  as  analyized  by  the 
Standards  Laboratory  associated  with  the  Chemical  Engineering 
Department  at  the  University  of  Alberta,  are  presented  on  Table  3.1. 
The  reported  properties  result  from  tests  performed  on  each  of  the 
oil,  prior  to  their  use  in  the  flume.  Extended  use  in  the 
experimental  study  resulted  in  slight  changes  in  the  physical 
properties  but  these  changes  were  ignored  in  the  following  analysis. 
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The  first  oil  studied  was  a  common  rapeseed  vegetable  oil 
available  from  any  retail  outlet.  This  oil  would  emulsify  rapidly 
during  the  tests  turning  the  water  into  a  milky  color  that  hindered 
the  observation  of  the  oil  slick.  In  addition,  this  oil  was 
unsuitable  for  the  laser  anemometer  velocity  measurements. 

The  second  oil  was  labeled  CWO  base  and  was  supplied  by  the 
Strathcona  refinery  of  Imperial  Oil  Company  Limited.  During  the 
testing  period  this  oil  remained  fairly  clear  and  as  a  result  it  was 
utilized  exclusively  in  the  velocity  measurements. 

The  last  oil  was  a  dark  heavy  machine  oil  acquired  from  dead 
storage  at  the  University  of  Alberta  Hydraulics  Laboratory.  Owing  to 
its  high  viscosity,  this  oil  was  utilized  in  the  investigation  of  the 
extreme  characteristics  of  viscosity. 

An  elaborate  oil  recovery  system  was  unwarranted  as  the 
amount  of  oil  (less  than  0.05  m^)  employed  during  individual  tests 
was  dimunitive.  Oil  was  introduced  into  the  flow  by  simply  pouring 
pre-measured  amounts  of  oil  upstream  of  the  slick.  The  oil  was 
allowed  to  flow  freely  from  a  bucket  onto  the  water  surface  in  an 
attempt  to  minimize  the  droplet  formation  and  entrainment. 

After  testing  was  completed,  the  oil  was  recovered  with  the 
use  of  an  ordinary  wet-dry  shop  vacumn  cleaner  and  stored  in  large 
plastic  containers  where  the  oil  and  water  mixture  from  the  vacumn 
cleaner  would  separate. 


CHAPTER  4 


ICE  COVERED  OIL  SLICKS 


4.1  General 


Previous  experimental  and  analytical  studies,  documented  in 
the  literature,  have  systematically  characterized  the  ice  free  oil 
slick  phenomena.  The  purpose  of  this  chapter  is  to  examine, 
analytically,  the  changes  induced  by  an  ice  cover  in  the  behaviour  of 
oil  slicks  retained  upstream  of  a  rigid  boom  in  flowing  water.  This 
treatment  will  be  limited  to  two  major  areas.  First,  the  addition  of 
a  ice  cover  will  alter  the  internal  velocity  distribution  of  the  oil. 
This  will,  subsequently,  affect  the  force  balance  of  the  oil  slick. 
The  second  major  effect  of  the  ice  cover  will  be  its  influence  on  the 
instability  of  the  oil-water  interface.  These  two  areas  will  be 
examined  in  detail  in  the  remainder  of  this  chapter. 

4.2  Velocity  Distribution 


Previous  investigators,  with  the  exception  of  Wicks,  have 
ignored  the  role  of  internal  oil  velocity  distributions.  The  momentum 
within  the  oil  phase  is  an  order  of  magnitude  less  than  the  underlying 
water  phase,  consequently  the  analysis  of  ice  free  oil  slicks  have  not 
been  adversely  hampered  by  neglecting  the  momentum  within  the  oil. 
However,  as  indicated  in  the  following  section,  the  oil  velocity 
distributions  within  the  ice  free  and  ice  covered  oil  slicks  are  shown 
to  be  markedly  different. 
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This  difference  stems  from  the  balance  of  mass  and  momentum 
in  the  oil  slick.  The  shear  force,  exerted  by  the  flowing  water,  on 
an  ice  free  oil  slick  induces  a  forward  (downstream)  flow  of  oil 
towards  the  boom.  However,  since  the  oil  is  arrested  at  the  boom,  a 
return  flow  of  oil  occurrs  in  the  upper  layer  of  the  slick.  The 
conservation  of  mass  dictates  that  this  mass  flux  of  the  return  flow 
must  equal  the  mass  flux  of  oil  in  the  downstream  direction.  The 
resulting  idealized  velocity  distribution  in  an  ice  free  oil  slick  is 
shown  in  Figure  4.1. a  as  a  typical  case. 

The  addition  of  an  ice  cover  over  an  oil  slick  restricts  the 
velocity  on  the  upper  surface  of  the  oil  to  zero,  in  order  to  satisfy 
the  condition  of  no  slip  at  a  (rigid)  boundary.  The  reduction  of  the 
velocity  profile  in  the  upper  layer  of  oil  to  zero  at  the  oil-ice 
interface  retards  the  return  momentum  flux.  However,  since  the  shear 
forces  on  the  oil-water  interface  induces  relatively  the  same 
downstream  flow  of  oil,  (a  similar  momentum  flux),  the  oil  slick  must 
thicken  to  compensate  for  the  momentum  deficit  resulting  from  the 
retarded  upper  layer  return  flow.  A  typical  idealized  velocity 
profile  in  an  ice  covered  oil  slick  is  shown  in  Figure  4 . 1 . b . 

The  idealized  velocity  distributions  can  be  developed  by 
solving  the  one  dimensional  Navier-Stokes  Equations.  It  will  be 
assumed  that  the  velocities  within  the  oil  are  laminar  in  magnitude 
and  the  flow  steady.  The  superimposed  orbital  velocities  induced  by 
the  passage  of  interfacial  waves  will  be  ignored. 
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Referring  to  the  idealized  velocity  distributions  depicted  in 
Figures  4.1. a  and  4 . 1 . b  ,  one  dimensional  flow  is  defined  when  the 
vertical  velocity  V  and  lateral  velocity  W  are  zero,  and  the 
longitudinal  velocity  U  is  a  function  of  longitudinal  distance  x  and 
vertical  distance  y  (ie.  U  =  u(x,y)).  In  addition,  the  pressure  P  is 
also  a  function  of  x  and  y  such  that  P  =  p(x,y).  These  conditions 
limit  the  general  form  of  the  Navier-Stokes  (N-S)  equations  to; 


N-S  Eqn.  x  direction 


-1  dP  .  d2U  d2U  1  dP 

=  T3xv'ay"  *  HF=u3x 


N-S  Eqn.  y  direction 


0  2 


dP 

3y 


N-S  eqn.  z  direction 


0  "37  >  p  -  P<x> 

Continuity  Eqn. 


(4.2.1) 


(4.2.2) 


=  0  »  U  =  u(y)  (4.2.3) 

Equations  (4.2.2)  and  (4.2.3)  indicate  that  the  piezometric 
pressure  is  a  function  of  x  only  and  velocity  is  solely  of  y, 
therefore  Equation  (4.2.1)  can  only  be  satisfied  if  the  left  hand  side 
and  the  right  hand  side  are  equal  to  a  constant  Cq. 


■ 
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1  dP  _  r 

u  (Jx  L0 


(4.2.4) 


The  general  solution  to  Equation  (4.2.4)  takes  the  following 

form; 


U  =  2^2  +  Cl*  +  C2 


(4.2.5) 


The  specific  solution  for  an  ice  covered  oil  slick  is  similar 
to  Couette  flow  between  two  parallel  flat  walls.  Referring  to  Figure 
4.1.b  the  following  boundary  conditions  are  evident. 


At  y  =  0  U  =  0  »  c2=0 


At  y  =  h  U  =  U-  interfacial  velocity 

■  I  _  J.  dP  h^_ 

h  2y  dx  h 

The  resulting  specific  solution  for  the  velocity  U  is; 


U 


y  U.  -n  dP  y  ,  ,  y  v 
=  ft  1  7T7  TT  V  1  “  J7  ) 


■dx  h 


(4.2.6) 


Schlichting  defines  a  dimensionless  pressure  gradiant  as; 


P*  =  h2  -dP  (4.2.7) 

2yUi  dx 

Substituting  the  dimensionless  pressure  gradiant  P*  -j n to 
equation  (4.2.6),  leads  to  the  following  dimensionless  velocity 
distribution  for  Couette  flow. 
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u  =  y  [  1  +  p*  (  1  +  y)] 


(4.2.8) 


An  additional  boundary  condition  can  be  derived  from 


continuity  of  mass.  It  was  stated  earlier  that  the  downstream  flow  of 
oil  equals  the  return  flow  of  oil  resulting  in  no  net  flow  of  oil  q 
past  a  given  section. 


(4.2.9) 


Performing  this  intergration  yields  a  dimensionless  pressure 
gradiant  of  P*  =  _3.  Thus,  Equation  (4.2.8)  becomes; 


(4.2.10) 


U  =  y  [  1  -  3  (  1  -  y) 
07  h  F 


IT 


The  maximum  upstream  velocity  of  U  =  -0.33  1).  occurs  at  y/h 
=  0.33  from  the  top.  At  y/h  =  0.67  the  velocity  is  zero. 


The  idealized  distribution  shown  in  Figure  4. 1 . b  represents 
Equation  (4.2.10). 

The  case  of  the  ice  free  oil  slick  velocity  distribution 
follows  the  preceeding  solution  but  with  a  different  boundary 
condition  on  the  free  surface.  Although  the  velocity  on  the  surface 
has  a  finite  value,  the  shear  stress  is  zero.  This  can  be  represented 
by  the  first  derivative  of  the  velocity  distribution  being  set  to  zero 


at  the  free  surface. 


' 
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At  y  =  0  dU  =  0  »  C1=  0 

cly  1 

Aty=h  y.=U  »  C9  =  U-  -  1  dP  h2 

1  2  1  27^ 

Utilizing  the  dimensionless  pressure  gradient  of  Equation 
(4.2.7)  the  specific  solution  of  Equation  (4.2.5)  in  a  dimensionless 
form  is; 


U  =  1  -  P*  (  1  -  y)  (4.2.11) 

F 

The  boundary  condition  of  zero  net  flow  past  a  given  section 
yields  a  dimensionless  pressure  gradiant  of  P*  _  _3/2.  Thus 
Equation  (4.2.11)  becomes; 

U  =1-3(1-  ,y,2)  (4.2.12) 

U7  1  W 

The  maximum  upstream  velocity  of  U  =  -0.5  now  occurs  at 
the  top  y/h  =  0  and  at  y/h  =  0.577  the  velocity  is  zero.  The 
idealized  velocity  distribution  presented  in  Figure  4.1. a  is  derived 
from  Equation  (4.2.12). 

4.3  Shear  Stress 


Development  of  interfacial  shear  stress  relationships  readily 
follows  the  preceding  laminar  flow  velocity  distribution  analysis. 
Laminar  shear  stress  is  proportional  to  the  first  derivative  of  the 
velocity  distribution  at  the  interface. 


q  >0  tm  ffct  tc  ***  nt‘MV)  K  ^yl 
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Differentiating  Equations  (4.2.10)  and  (4.2.12)  yields  the 
following  shear  stress  relationships; 

Ice  free  top  y  =  0  ot  =  0 

interface  y  =  h  <j.=  y  3U^  •  (4.3.1) 

~TT 

Ice  covered  top  y  =  0  ot=  -y  2U^  (4.3.2) 

“FT 

interface  y  =  h  a.  =  y  4U^  (4.3.3) 

F~ 

A  comparison  can  be  made  concerning  the  relative  oil  slick 
thickness  between  the  ice  free  and  ice  covered  conditions.  Assuming 
that  for  similar  approach  velocities  the  interfacial  velocities  at  a 
given  distance  from  the  leading  edge  of  the  slick  will  be  basically 
the  same  for  ice  free  and  ice  covered  oil  slicks,  the  interfacial 
shear  stresses  should  be  similar.  Equating  Equations  (4.3.2)  and 
(4.3.3)  with  similar  interfacial  velocities  and  shear  stress  results 
in  the  thickness  of  an  ice  covered  oil  slick  4/3  times  the  thickness 
of  an  ice  free  oil  slick. 

The  utilization  of  shear  stress  estimates  can  be  employed  in 
the  development  of  force  balance  relationships  of  the  oil  slick. 

Consider  the  ice  free  oil  slick  control  volume  presented  in 
Figure  4.2. a,  the  sum  of  horizontal  force  states; 

F 1  +  P  3  +  a-j  AX  = 


(4.3.4) 
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or 


lp0gh2  +pQgpw(h+Ah)Ah  +o.ax  =l_pQgh2  +pQg(h+Ah)Ah  (4.3.5) 


Assuming  Ah/2  <<  h  and  rearranging  Equation  (4.3.5)  yields; 


O.  AX  =  pQ  g  (1  -  P0)  h  Ah  (4. 3. 6. a) 

In  terms  of  the  interfacial  slope  Am  this  is 


a i  ax  =  pw  g  (1  -  pQ)  h  Am  (4.3.6.b) 

Referring  to  Figure  4 . 2 . b  the  force  balance  for  a  control 
volume  of  oil  under  an  ice  covered  oil  slick  leads  to  a  similar 
solution  in  which; 


(a  -j  -  °  t )  AX 


pw  g  (i  -  p0)  h  Ah 


w 


(4.3.7) 


Equation  (4.3.7)  applies  equally  to  the  increase  in  oil 
thickness  and  the  interfacial  slope. 


Utilization  of  the  previously  developed  shear  stress 
relationships  introduces  a  unique  relationship. 

For  the  ice  free  case,  substituting  Equation  (4.3.1)  for  the 
interfacial  shear  stress  of  Equation  (4.3.6)  yields; 
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3  y  U-  AX 
TT“ 


p0  g  (1  -  pQ)  h  ax 


or  rearanging 

-3  =  -h2  [Pq  g  (1  -  pQ)  Ah 

2  2  UU.  ~  AX 

l  w 


(4.3.9) 


(4.3.10) 


This  corresponds  directly  to  Schlichting  dimensionless 
pressure  gradiant  equation  (4.2.7)  in  which. 


dP  =  p0  9  U  -  p0)  Ah  (4.3.11) 

(lx  -  AX 

PW 

and  results  in  the  same  Dimensionless  pressure  gradiant 
P  =  -3/2  as  given  by  velocity  solution  of  N.S.  Equations. 

In  a  similiar  manner  substitution  of  Equation  (4.3.2)  and 
(4.3.3)  for  the  interfacial  and  top  shear  stress  into  Equation  (4.3.8) 
yields; 

(4.3.12) 

(4.3.13) 


As  earlier  presented  the  Dimensionless  pressure  gradiant  is 
now  P  =  -3. 


-6  ii  Ui  =  pw  g  (1  -  pQ)  h  Ah 


IT 


w 


or 

-3  = 


-  h  p  g  (1  -  p  )  Ah 

2  y  U-  0  —  AX 

l  p 

w 


•  •  • 
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The  direct  consequence  of  these  pressure  gradiants  can  be 
shown  by  the  comparison  of  the  interfacial  slopes  on  ice  covered  and 
ice  free  oil  slicks.  The  two  slopes  are  related  as; 

Ah  ice  covered  =  2.0  p  All  ice  free  (4.3.14) 

ax  —  AX 

pw 

In  conclusion,  the  difference  between  the  velocity 
distributions  and  related  force  balance  of  ice  covered  and  ice  free 
oil  slicks  plays  a  dominant  role  in  describing  the  influence  of  an  ice 
cover  on  an  oil  slick 

4.4  Interfacial  Instability 


It  has  been  noted  by  various  researchers  that  the  limiting 
factor  to  successful  containment  of  oil  slicks  is  the  onset  of 
breaking  interfacial  waves  at  water  velocities  of  the  order  of  0.4  m/s. 
Section  2.4  on  interfacial  instabilities  indicates  that  the 
development  and  eventual  breakup  of  interfacial  wave  forms  can  be 
described  by  Kelvin  -  Helmholtz  (K-H)  instability  analysis.  Even 
though  the  solutions  to  linear  K-H  analysis  do  not  predict  the  exact 
magnitude  of  wave  breakup  K-H  analysis  is  sufficiently  capable  of 
indicating  trends  and  order  of  magnitude  on  interfacial  wave 
instabi 1 i ti ty. 

One  of  the  earlier  stated  objectives  was  to  document  the 
development  of  interfacial  waves  as  influenced  by  an  ice  cover.  The 
following  section  outlines  the  development  of  K-H  instability  analysis 
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as  governed  by  the  addition  of  a  solid  cover. 

Assuming  the  flow  in  the  water  and  oil  phases  is 
irrotational ,  the  velocity  potentials  and  small  perturbations  must 
satisfy  the  Lapace  Equations.  Denoting  the  interface  between  the  oil 
and  water  as  y  =  0,  the  following  two  boundary  conditions  (B.C.)  exist; 


B.C.  1)  At  the  ice  cover,  y  =  h,  the  velocity  normal  to  the 

ice  is  zero.  Thus  VQ  =  0  =  d^ 

dy” 

B.C.  2)  If  the  bottom  of  the  water  is  suficiently  far,  no 
perturbations  exist.  Thus  0  •  _  o  at  y  =  -infinity 

The  general  solutions  to  the  Lapace  Equations  with  the 
preceding  boundary  conditions  are; 


0Q'  =  A0  cosh  K(y-h)  exp  [i(kx-at)] 

V  =  Aw  exp  Ky  exp  [ i (kx-at )] 

In  which, AQ  s  aw  are  constants 
K  =  Wave  number  =  2W1 
a  =  wave  period  =  2n-/1 


(4.4.1) 


(4.4.2) 


The  interfacial  wave  form  of  amplitude  a,  can  be  described  as; 
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n  =  a  exp  [i(Kx-ot)]  (4.4.3) 

The  total  derivative 


D_  (y-n(x.t) )  =  0  (4.4.4) 

Dt 

states  that  a  particle  on  the  interface  remains  on  the  interfacial 
streamline.  Expanding  Equation  (4.4.4)  and  evaluating  at  the 
interface  y  =  0  (Taylors  expansion)  yields  the  following  two  kinematic 
boundary  conditions; 


B.C.  3)  dn  +  U  dn 
dt  0  3^ 


B.C.  3)  dn  +  U  dn 
<Jt  w  3x 


4' 

12. 

dy 

4' 
_ w 

dy 


at  y  =  0 


at  y  =  0 


(4.4.5) 


(4.4.6) 


An  additional  dynamic  boundary  condition  is  derived  from  the 
balance  of  pressures  at  the  interface  with  surface  tension  such  that; 


B.C.  5)  Pw  -  PQ  =  -T  d2n  at  y  =  0 

d? 

The  interface  is  a  common  streamline  between  the  two  fluids, 
and  thus  the  pressure  can  be  evaluated  with  a  linearized  form  of  the 
Bernoulli  Equation  and  Boundary  Condition  5  becomes; 


pw(d^  +9n)  -  Po(d0;  +UQd^  +gn)  =  -T  d_n 

dt  dx  dt  dx  dx^ 


(4.4.7) 
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Substituting  the  general  solution  to  the  Lapace  Equations 
(4.4.1)  and  (4.4.2)  and  the  assumed  interfacial  shape.  Equation 
(4.4.3)  into  the  Boundary  Conditions  3,  4,  and  5  (Equations  (4.4.5), 
(4.4.6),  and  (4.4.7))  results  in  the  following  three  equations; 

3  i(U0K  -  a)  =  K  A0  sinh  ( Kh )  (4.4.8) 

a  i(UwK  -  a)  =  _K  (4.4.9) 

°wCAwi(uwK-<J)+ga]  -P0[A0i  cosh  Kh  (U0K-o)+ga]  =  TaK2  (4.4.10) 

The  simultaneous  solution  for  the  three  unkowns  a,  A  and 
Aw  yields  the  following 


pw(uw-c)2  +p0(u0-c)2  coth  Kh  ~TK  -a  (h„-p0)  =  o  (4.4.11) 

K 


Solving  for  the  interfacial  wave  speed  C  yields 


C  =  p  U  +U  D  ±  [g(p  -p  )+TK2  -p  D  (U  -U  )2]  1/2 
Mw  w  o  LavMw  woJ  Mw  v  w  o'  J 


(4.4.12) 


p  +  D 
Mw 


K  pw  +  D 


(»w+  °r 


In  which  D  =  pQcoth  (Kh) 


The  first  term  of  Equation  (4.4.12)  is  the  weighted  mean 
velocity  of  flow  and  the  terms  within  the  square  root  indicate  the 
interfacial  wave  velocities.  Instabil itity  of  the  interface  will 
occur  when  C  becomes  imaginary.  Denoting  the  critical  flow  velocity 


-  Uw  -  U 


0 


» 


instability  is  defined  as; 
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Ucr  =  (»«)  po  coth  Kh  [  a  (ow  -  °o)  +  Tl<]  (4.4.13) 

pw  +  P0  coth  Kh  K 

The  density  of  water  in  the  S.I.  system  is  1  g/cnP,  thus 
for  convenience  Equation  (4.4.13)  reduces  to; 

uCp  =  1  +  tanh  Kh  [  g  (1  -  pq)  +  TK  ]  (4.4.14) 

po  K 

It  is  evident  from  Equation  (4.4.14)  that  for  a  given  wave 
Number  K  the  critical  velocity  Ucr  increases  with  the  depth  of  oil  h. 

Reduced  to  a  similar  form,  the  analysis  of  ice  free 
interfacial  waves  presented  by  Leibovich  (Section  2.4  and  Equations 
(2.4.40)  and  (2.4.42))  becomes; 

For  thin  si icks  as  h  >  0 

Ucr  =  9/K  +  K  (T0/a  +  To/w)  (4.4.15) 

and  for  thick  slicks  as  h  >  infinity 

Ucr  =  U  +  Up0)[g/K  (1  -  p0)  +  KT0/W]  (4.4.16) 

* 

As  previously  stated  Equations  (4.4.15)  and  (4.4.16)  show 
that  ice  free  oil  slicks  are  destabilized  by  increasing  thickness. 

Inti  ally  it  would  appear  that  since  Equation  (4.4.14)  indicates  the 
converse  for  an  ice  covered  oil  slick  that  this  alone  should  suffice 
in  describing  the  stability  of  ice  covered  oil  slicks.  However,  close 
examination  reveals  that  for  thick  slicks.  Equation  (4.4.14)  equals 
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Equation  (4.4.16)  and  both  ice  free  and  ice  covered  thick  slicks  have 
similar  instability.  The  conclusion  drawn  from  a  comparison  of  all 
three  equations  would  be  that  ice  free  slicks  should  be  as  stable,  if 
not  more,  than  ice  covered  oil  slicks. 


CHAPTER  5 


EXPERIMENTAL  RESULTS 


5.1  General 


The  intent  of  this  chapter  is  to  present  the  results  of  the 
experimental  phase  of  the  study.  As  indicated  earlier,  the  objectives 
of  this  investigation,  aside  from  the  initial  familiarization  with  the 
behaviour  of  an  ice  free  oil  slick,  were  three-fold.  First,  a  series 
of  tests  was  conducted  in  conjunction  with  the  control  runs  to 
identify  the  basic  behavioural  differences  resulting  from  the  addition 
of  an  ice  cover.  Secondly,  the  distribution  of  internal  oil 
velocities  was  investigated  both  with  and  without  an  ice  cover. 
Finally,  the  formation  of  interfacial  waves  as  influenced  by  the 
addition  of  an  ice  cover  was  examined. 

5.2  Test  Schedule 


The  afore-mentioned  objectives  were  realized  through  three 
series  of  tests.  The  three  series  are  broadly  defined  in  the 
following  paragraphs.  It  should  be  noted  at  this  point  that 
observations,  documentation  and  objectives  of  several  tests  overlapped 
the  series  boundaries. 
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Series  A  Preliminary  Testing. 

This  series  served  to  produce  control  tests  of  ice  free  oil 
slick  behaviour  and  to  develop  an  intial  understanding  of  the 
influence  of  a  solid  ice  cover.  Series  A  was  conducted 
exclusively  with  the  rapeseed  oil.  A  total  of  15  tests 
consisted  of  longitudinal  profile  measurements  were  performed. 

Series  B  Velocity  Distributions. 

This  series  was  restricted  to  defining  the  internal  oil  slick 
velocity  distributions  in  both  ice  covered  and  ice  free  oil 
slicks.  A  total  of  9  runs  were  conducted  of  which  8  utilized 
the  CWO  base  oil,  and  1  used  rapeseed  oil.  In  addition  to 
the  velocity  data,  longitudinal  profiles  were  recorded. 

Series  C  Interfacial  Waves. 

This  series  of  testing  consisted  of  filming  20  individual 
tests.  Stop  action  analysis  of  the  films  facilitated 
documentation  of  interfacial  waves,  lgngitudial  profiles,  and 
observation  of  the  headwave  developments.  Series  C  involved 
testing  of  all  three  oils. 
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Table  5.1  summarizes  the  complete  range  of  test  conditions  utilized  in 
the  44  tests. 

The  synopsis  of  test  results  are  subsequently  presented  as  a 
general  discussion  preceding  the  individual  interpretation. 

5.3  Series  A  Test  Results 


A  brief  description  of  the  overall  oil  slick  behaviour  with 
and  without  an  ice  cover  is  a  preliminary  requisite  for  the 
appreciation  of  the  experimental  results.  Information  contained  in 
this  section  is  primarily  a  capsule  summary  of  the  observations 
concluded  during  the  course  of  Series  A  testing.  However,  some  of  the 
general  observations  pertaining  to  different  oils  (ie. Series  B  and  C) 
are  included  in  this  section.  Subsequent  sections  will  present  the 
detailed  experimental  results  of  Series  B  and  C. 

Basically,  a  rigid  boom  of  sufficient  draft  was  deployed  into 
the  flume  to  arrest  and  pond  a  fixed  volume  of  oil  in  order  to  form  an 
oil  slick  of  some  deterministic  shape.  The  longitudinal  profiles 
corresponding  to  the  44  tests  are  recorded  in  Appendix  A.  The 
resulting  shape,  the  oil  slick  assumed,  can  be  characterized  by  three 
velocity  ranges.  The  transition  of  the  shape  between  these  ranges  was 
not  clearly  distinguishable  and  the  range  limits  vary  with  the 
properties  of  the  three  oils  tested. 
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Typically,  at  low  water  velocities  of  the  order  of  less  than 
0.2  m/s,  the  oil  slick  assumed  a  smooth  parabolic  shape  with  little  or 
no  evidence  of  either  a  headwave  or  interfacial  waves. 

As  the  water  velocities  were  increased  to  somewhere  between 
0.2  and  0.4  m/s,  the  oil  slick  profile  obtained  a  more  distinct  shape, 
characterized  by  the  formation  of  a  headwave  at  the  leading  edge  of 
the  slick  and  the  development  of  interfacial  waves  along  the  entire 
length  of  slick.  For  a  given  volume  of  oil,  the  oil  slick  length 
shortened  as  the  velocity  was  incremented. 

Finally,  the  last  stage  was  obtained  as  the  water  velocities 
were  further  increased  to  exceed  approximately  0.4  m/s.  The  oil  slick 
continued  to  thicken  and  decrease  in  length  until  the  interfacial 
waves  started  to  break  in  the  vicinity  of  the  headwave.  Oil  droplets 
were  entrained  into  the  flow  and  when  water  velocities  were 
sufficiently  high,  the  oil  droplets  did  not  rise  back  and  coalesce 
with  the  slick  but  rather  swept  past  the  boom.  The  breakup  of 
interfacial  waves  resulted  in  a  rapid  loss  of  oil  past  the  boom. 
Consequently,  oil  profiles  and  measurements  in  this  study  were 
restricted  to  velocities  below  this  range. 

The  preceding  description  applies  equally  to  the  ice  free  and 
ice  covered  oil  slicks,  although  the  range  limits  varies  between  the 
three  distinguishable  stages. 
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It  was  observed  that  ice  covered  oil  slicks  exhibit  two  major 
deviations  from  their  ice  free  counter-parts.  The  most  obvious 
difference  is  that  for  similiar  relative  velocities,  ice  covered  oil 
slicks  are  thicker  and  shorter  than  ice  free  oil  slicks.  In  addition 
the  headwave  itself  appears  to  be  thicker  but  less  pronouced  at  the 
downstream  contraction  than  at  similar  velocities  in  the  ice  free 
case.  The  second  dissimilarity  is  the  apparent  stability  of  the 
interfacial  waves  and  the  associated  diminished  entrainment  failure  at 
comparable  approach  velocities  to  those  of  the  ice  free  oil  slicks. 

As  indicated  above,  ice  covered  oil  slicks  are  shorter  than 
their  ice  free  counter-parts  at  similar  approach  velocities.  The  ice 
covered  oil  slick  exhibit  an  analogous  behaviour  to  ice  free  oil  slick 
in  their  tendency  to  shorten  with  increasing  approach  velocity. 

In  order  to  represent  the  various  oil  slick  lengths  as  a 
function  of  approach  velocity  a  normalized  "effective  length"  is 
defined.  The  effective  length  directly  incorporates  the  relative 
influence  of  headwave  thickness,  near  boom  region,  volume  of  oil,  and 
the  flume  width. 

* 

The  first  factor  to  be  considered  in  defining  an  effective 
length  is  the  intial  thickness  of  the  headwave.  The  relative  headwave 
thickness  of  two  different  oils  should  be  related  by  the  ratio  of 
their  specific  weights  (Equation  2.2.1).  Drawing  upon  the  difference 
between  Benjamin  and  Wilkinson's  definition  of  the  Headwave  Froude 
Number,  the  ice  covered  headwaves  exceeds  the  thickness  of  ice  free 
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headwaves  by  the  inverse  of  specific  weights.  The  definition  of 
effective  Length  uses  the  specific  gravity  of  rapeseed  oil  to 
normalize  all  44  tests. 

A  second  factor  in  the  definition  of  the  effective  length 
pertains  to  the  influence  of  the  near  boom  region.  In  the  immediate 
vicinity  of  the  boom  the  thickness  of  oil  is  primarly  governed  by  the 
stagnation  pressure  induced  by  the  protruding  boom  draft.  In  this 
case  it  was  assumed  that  this  influence  occurs  for  a  reach  length 
upstream  of  the  boom  equal  to  5.0  times  the  boom  draft.  Consequently 
5.0  boom  drafts  are  subtracted  from  the  measured  oil  slick  length. 

Finally  the  length  of  an  oil  slick  is  governed  by  the  volume 
of  oil  tested  and  the  width  of  the  flume.  In  this  case  the  width  of 
the  flume  is  constant  for  all  tests  and  is  not  included  in  the 
effective  length. 

Coupling  the  preceeding  factors  together,  the  effective 
length  of  the  oil  slick  is  defined  as; 

Ice  free; 
s 

El.  =  rape  length  -  5.0  boom  draft 
Soi  1 

volume 

Ice  covered; 

s  , 

El.  =  rap-e-  x  - x  length  -  5.0  boom  draft 

soi 1  so i 1 _ 

volume 
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Utilizing  these  relationships,  all  44  test  have  been  plotted 
as  Effective  length  (El.)  vs  Approach  Velocity  and  are  presented  in 
Figure  5.1  . 


5.4  Series  B  Test  Results 


Table  5.1  reveals  that  9  runs  were  conducted  to  investigate 
velocity  distributions.  A  total  of  10  velocity  profiles  were  recorded 
of  which  4  were  under  an  ice  cover.  These  10  velocity  profiles  are 
plotted  in  Appendix  B 

In  order  to  validate  the  theoretical  velocity  distributions 
developed  in  Chapter  4  , Equation  (4.2.10)  (ice  covered)  and  Equation 
(4.2.12)  (ice  free),  it  is  neccesary  to  reduce  these  profiles  into  a 
non-dimensional  form. 

As  evident  in  Appendix  B  the  interfacial  velocities  were  not 
directly  measured  by  the  laser  anemometry  system.  The  interface 
oscillated  verticaly  with  each  passage  of  an  interfacial  wave  and 
since  the  laser  was  aimed  at  a  fixed  vertical  location,  the  readings 
were  conducted  in  a  continuously  changing  medium  between  water  and 
oil.  The  coefficients  utilized  with  the  laser  system  in  calculating 
velocity  differ  with  oil  and  water,  and,  thus,  any  attempt  to  directly 
calculate  the  interfacial  velocity  was  impractical. 
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Estimates  of  the  interfacial  velocities  were  derived  from  the 
internally  measured  velocities  and  with  the  aid  of  the  preceding 
theoretical  relationships.  In  the  case  of  ice  free  conditions,  the 
surface  velocity  was  first  established  by  extrapolating  the  measured 
velocities  immediately  below  the  surface.  This  velocity  was  in  turn 
assumed  to  be  -0.5  times  the  interfacial  velocity  as  suggested  by 
Equation  (4.4.12).  For  the  ice  covered  condition,  the  interfacial 
velocity  was  calculated  as  -3.0  times  the  maximum  upstream  velocity. 
Both  the  surface  and  interfacial  velocities  developed  in  this  manner 
are  indicated  in  Appendix  B. 

Utilizing  the  estimates  of  interfacial  velocity, 
dimensionless  velocity  profiles  composed  of 

U(measured)/U.|  (calculated)  vs  y/h  were  established.  The  profile  for 
the  ice  free  and  ice  covered  oil  slicks  are  presented  in  Figures  5.2 
and  5.3  respecti vely. 

The  development  of  the  internal  velocity  distributions  in 
Chapter  4  leads  to  an  examination  of  shear  stress  and  a  comparison  of 
the  relative  slick  thickness  between  ice  covered  and  ice  free  oil 
slicks.  It  was  demonstrated  that  for  similar  shear  stresses  and 
interfacial  velocities  the  ice  covered  oil  slick  should  be  4/3  times 
as  thick  as  their  ice  free  counterpart.  Visual  comparison  of  the  oil 
profiles  contained  in  Appendix  A  for  similar  approach  velocities 
reveal  that  this  trend  is  reasonable.  Specific  comparison  is  not 
possible  as  the  thickness  of  oil  slicks  varies  with  the  longitudinal 
position  and  the  length  of  ice  free,  in  addition,  ice  covered  oil 
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slicks  also  vary. 

Support  of  the  assumption  of  a  similar  interfacial  velocity 
and  shear  stress  is  provided  by  the  following  analysis.  Figure  5.4 
represents  the  interfacial  velocity  as  a  function  of  approach  velocity 
for  both  ice  free  and  ice  covered  oil  slicks.  The  general  trend 
indicates  both  ice  free  and  ice  covered  interfacial  velocities  are 
similar.  Figure  5.5  reveals  the  shear  stress  Vs.  U-j;  as 
calculated  by  Equation  (4. 3.1), for  the  ice  free  case  ,  and  (4.3.3)  for 
the  ice  covered  example.  Here  again,  the  ice  free  and  ice  covered 

represent  similar  shear  stress  at  comparable  velcoities.  The 
preceding  two  plots  do  not  distinguish  as  to  the  longitudinal  location 
of  the  recorded  velocity  data.  However,  Figure  5.6  reveals  the  shear 
stress  coefficient  as  a  funciton  of  longitidinal  Reynolds  Number.  The 
coefficient  is  calculated  as; 


2  a 
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Table  5.2  reveals  the  primary  data  used  to  compile  the 
preceding  Figures  5.4,  5.5,  and  5.6  . 


5.5  Series  C  Test  Results 

It  was  observed  in  ice  free  oil  slicks  that  intially  short 
waves,  at  the  leading  edge,  of  the  order  of  5  cm.  would  grow  as  they 
migrated  along  the  slick  to  become  larger  waves  of  the  order  15  -  30 
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cm.  within  the  viscous  zone.  These  waves  would  travel  downstream  with 
the  following  speed; 

Rapeseed  oil  15  -  25  cm/s. 

CWO  base  oil  10  -  20  cm/s. 

Heavy  gear  oil  3-6  cm/s. 

Since  a  complete  range  of  superimposed  waves  existed  along  the  slick, 
encompassing  various  wave  lengths  and  periods,  it  was  extremely 
difficult  to  pin  point  specific  characteri sties  or  develop 
relationships  for  the  waves  as  functions  of  either  the  different  oils 
or  approach  velocities  in  the  various  runs. 

However,  the  major  observation  noted,  though  not  quantified, 
was  the  addition  of  an  ice  cover  markedly  reduced  the  interfacial  wave 
activity.  Interfacial  waves  under  an  ice  cover  tend  to  travel  slower 
and  are  composed  of  only  a  few  waves  with  small  amplitudes.  It  was 
noted  that  the  waves  consisted  of  a  longer  wave  length  than  those  in 
the  ice  free  oil  slick.  The  overall  basic  observation  during  the 
study  indicated  that  ice  covered  oil  slicks  exhibit  decreased  wave 
activity,  a  stable  interface,  and  are  less  susceptible  to  entrainment 
failure  at  comparable  velocities  to  ice  free  oil  slicks. 
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CHAPTER  6 


DISCUSSION 


6.1  General 


The  observations  of  ice  free  oil  slicks  made  during  this 
study  did  not  deviate  to  any  significant  degree  from  those  reported  in 
the  literature  reviewed.  This  chapter  primarily  addresses  the 
observed  behaviour  of  oil  slicks  as  influenced  by  an  ice  cover. 

The  general  characteristics  of  ice  covered  oil  slicks 
simulate  those  of  ice  free  oil  slicks.  For  example,  both  types  of  oil 
slicks  exhibit  a  decrease  in  oil  slick  length  with  an  increase  in 
velocity  as  indicated  in  the  effective  length  plot  of  figure  5.1  . 

With  the  exception  of  the  heavy  oil  (in  which  the  viscosity  greatly 
exceeds  the  other  oils),  this  type  of  dimensional  length/volume 
adequately  describes  the  general  length  relationship  as  a  function  of 
approach  velocity  for  both  ice  free  and  ice  covered  oil  slicks. 
Deviations  from  a  single  line  arise  from  experimental  error,  viscosity 
effects,  and  error  in  definition  of  an  effective  length  (ie.  boom 
draft  factor). 

The  three  major  differences  in  behaviour  of  ice  covered  oil 
slicks  are  the  altered  headwave  activity,  a  general  increased  oil 
slick  thickness,  and  the  apparent  stability  of  the  oil-water 
interface.  These  variations  will  compose  the  leading  theme  of  the 
following  sections. 
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6.2  Headwave 


The  development  of  the  headwave  is  the  one  of  the  most 
dominant  features  of  oil  slicks  retained  by  a  boom.  A  simple  Headwave 
Froude  Number  states  that  the  thickness  of  the  headwave  is  twice  the 
thickness  of  the  slick  in  the  lee  of  the  headwave.  More  specifically, 
the  thickness  is  a  function  of  Froude  Number  and  relative  depth. 
Benjamin  argued  that  for  conservation  of  energy  and  momentum  the 
formation  of  the  headwave  is  imperative.  On  the  other  hand,  Milgram 
and  Van  Houten  contradicted  this  premise  by  arguing  that  the  shear 
stresses  are  all  important  in  the  formation  of  the  oil  slick  profile. 
They  further  demonstrate  (at  low  velocities)  that  in  the  absence  of 
dynamic  considerations  the  shear  stress  can  adequately  describe  the 
oil  slick  shape.  However,  at  higher  velocities,  when  the  headwave  is 
visually  obvious  its  existence  can  not  be  described  from  shear 
stresses  alone.  It  was  observed  during  the  course  of  this  study  that 
with  increasing  velocity  the  headwave  was  more  evident. 

It  is  reasonable  to  consider  that  both  of  these  arguments  are 
valid  to  a  degree.  At  lower  velocities,  the  shear  stresses  are 
capable  of  consuming  the  energy  necessary  to  satisfy  Benjamin's 
momentum  requirements,  thus  explaining  the  absence  of  any  headwaves  in 
the  low  velocity  ranges.  At  velocities  exceeding  about  0.2  m/s,  the 
headwave  must  form  to  meet  the  deficiency  in  energy  consumption  not 
supplied  by  shear  stresses.  Observations  within  this  study  indicate 
that  the  headwave  only  aproached  2  times  the  thickness  at  the  higher 
end  of  the  velocities  tested. 
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It  was,  also,  observed  that  the  ice  covered  oil  slick 
exhibits  a  thicker  headwave  but  it  is  less  pronounced  than  in  the  ice 
free  state  at  the  downstream  contraction.  This  could  be  caused  by  two 
factors.  First,  the  velocity  field  approaching  the  oil  slick  is 
somewhat  reduced  in  the  boundary  layer  induced  by  the  presence  of  an 
ice  cover,  possibly  explaining  the  less  pronounced  headwave  formation. 

The  second  possible  factor  describing  the  increased  thickness 
can  be  explained  by  comparing  the  analysis  of  headwaves  by  Benjamin 
and  Wilkinson.  Wilkinson  allows  the  oil  slick  to  rise  above  the  water 
and  assumes  the  stagnation  point  at  the  water  surface.  However, 
Benjamin's  analysis  of  air  cavities  assumes  the  stagnation  point  is 
adjacent  the  cavity  roof.  This  corresponds  directly  to  an  oil  slick 
under  an  ice  cover.  Close  examination  of  the  two  headwave  thickness 
relationships  reveals  that  Benjamin's  thickness  exceeds  that  of 
Wilkinson  by  the  factor  1/s  (this  factor  was  utilized  in  the  defintion 
of  the  effective  length).  Thus,  as  observed,  it  is  valid  that  the 
headwave  is  thicker  under  an  ice  cover. 

6.3  Velocity  Distribution 


The  dimensionless  velocity  profiles  presented  in  Chapter  5 
(Figures  5.2  and  5.3)  appear  to  be  in  reasonable  agreement  with  the 
analytical  velocity  distributions  developed  in  Chapter  4. 

Examination  of  the  velocity  profiles  contained  in  Appendix  B 
reveal  that  the  magnitude  of  velocity  is  such  that  the  flow  is 
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basically  laminar.  This  supports  the  use  of  the  simple  1-d  solution 
of  the  Navier  -  Stokes  Equations  for  the  velocity  distributions. 

Visual  comparison  of  the  two  distributions  suggests  that  the 
deviations  from  the  theoretical  relationship  are  greater  for  the  ice 
free  case.  This  deviation  from  the  theory  probably  stems  from  two 
sources.  First,  since  the  interfacial  waves  migrate  downstream  and 
induce  the  oil  slick  to  oscillate  vertically,  the  fixed  position  laser 
anemometry  measurements  are  actually  averaging  the  velocities  from 
several  different  oil  depths.  The  interfacial  waves  were  most  active 
in  the  ice  free  oil  slick  since  the  addition  of  an  ice  cover 
suppressed  wave  activity.  The  second  source  of  error  is  from  the 
neglect  of  the  orbital  velocities  induced  by  the  passage  of 
interfacial  waves.  This  error  would  be  significant  near  the  interface 
and  of  even  more  importance  in  the  ice  free  state  with  its  increased 
wave  activity.  Evidence  of  this  error  is  indicated  by  the  slower 
interfacial  velocities  measured  by  floating  particles  as  shown  in 
Appendix  B.  The  particle  velocities  at  the  interface  were 
considerably  slower  than  the  measured  surface  velocities  and  therefore 
demonstrated  the  trend  for  greater  error  close  to  the  interface.  In 
summary,  even  though  some  deviations  from  theory  is  evident,  the 
general  agreement  between  the  measured  and  theoretical  velocites 
appears  to  be  acceptable. 

In  conclusion,  the  difference  between  the  velocity 


distributions  of  ice  covered  and  ice  free  oil  slicks  plays  a  dominate 
role  in  describing  the  influence  of  an  ice  cover  on  an  oil  slick. 
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6.4  Interfacial  Instability 

The  second  major  deviation  in  the  oil  slick  behaviour, 
induced  by  the  presence  of  an  ice  cover,  is  the  apparent  suppression 
of  interfacial  wave-  growth. 

There  are  a  several  possible  reasons  as  to  why  the  preceding 
conclusions  are  in  direct  contradiction  to  the  linear  K-H  instablitity 
analysis.  First,  it  has  been  suggested  that  interfacial  wave  growth 
is  suppressed  by  the  presence  of  a  thick  boundary  layer.  In  the  case 
of  an  ice  covered  oil  slick  the  boundary  layer  already  present  from 
the  ice  cover  itself  should  suppress  interfacial  waves  at  the 
interface.  This  is  supported  by  the  observation  of  waves  of  greater 
length  on  an  ice  covered  oil  slick  which  are  less  critical  than  the 
shorter  waves  that  lead  to  breakup  in  K-H  instability  analysis. 

Secondly,  viscosity  dampens  the  growth  of  interfacial  waves. 
The  increased  thickness  of  the  oil  slick  as  previously  demonstrated 
allows  a  greater  internal  dissipation  of  kinetic  energy.  This 
mechanism  for  viscous  energy  dissipation  appears  to  suppress  the  K-H 
instability  at  the  oil-water  interface. 

Finally,  the  preceding  analysis  has  been  simplified  in 
several  other  areas.  The  K-H  solution  is  based  on  the  assumption  of 
irrotational  flow.  The  analytical  and  measured  velocity  distributions 
indicate  that  this  is  not  valid  for  either  the  ice  free  and  ice 
covered  cases.  Liebovich  shown  that  neglect  of  the  non-linear  terms 
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in  the  solution  of  K-H  instability  analysis  lead  to  lower  critical 
velocities  than  those  which  exist  before  the  onset  of  interfacial  wave 
break-up. 


In  conclusion,  although  K-H  instablity  analysis  indicates 
that  ice  free  oil  slicks  should  be  stable  in  comparison  to  ice  covered 
oil  slicks,  the  observations  during  this  study  support  the  preceding 
arguments  that  ice  covered  oil  slicks  are  less  subject  to  instability 
failure  than  their  ice  free  counterparts. 


CHAPTER  7 


CONCLUSIONS  AND  RECOMMENDATIONS 


7.1  CONCLUSIONS 


The  general  behaviour  of  ice  free  oil  slicks  is  well 
documented  in  the  available  literature.  This  study  serves  as  an 
initial  evaluation  of  the  effects  a  solid  ice  cover  causes  on  the 
development  of  an  oil  slick  retained  by  a  boom.  The  major  conclusions 
arising  from  this  study  are  as  follows: 

1.  There  were  no  deviations  in  the  behaviour  of  ice  free  oil 
slicks  which  were  contary  to  the  literature. 

2.  The  headwave  under  an  ice  cover  was  observed  to  be 
thicker  than  those  of  ice  free  oil  slicks 

3.  Unique  relationships  for  the  internal  oil  velocities  were 
introduced.  The  differences  between  the  ice  covered  and  ice 
free  oil  slick  velocities  play  a  dominate  role  in  describing 
the  influence  of  an  ice  cover. 

4.  As  a  result  of  the  differences  in  velocity  disrtibutions, 
the  ice  covered  oil  slicks  are  thicker  and  shorter. 
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5.  A  linear  Kelvin  -  Helmholtz  instability  analysis  indicate 
the  stability  of  ice  covered  oil  slicks  increases  with 
increasing  oil  slick  thickness.  However,  K-H  instability 
analysis  shows  that  ice  free  oil  slicks  should  always  be  more 
stable  than  ice  covered  oil  slicks.  The  observations 
conducted  in  this  study  lead  to  a  direct  contradition  to  the 
K-H  instablity  analysis  as  ice  covered  oil  slicks  are  more 
stable  than  their  ice  free  counterparts. 

7.2  RECOMMENDATIONS 


Prior  to  the  establishment  of  a  comprehensive  analytical 
model  of  ice  covered  oil  slicks,  several  critical  areas  should  be 
examined; 


1.  In  order  to  facilitate  the  use  of  shear  stress  in  a  force 
balance,  a  more  detailed  velocity  survey  should  be  conducted 
along  the  entire  oil  slick  length. 

2.  The  question  of  interfacial  instablity  shoud  be 
addressed.  Emphasis  on  the  establishment  of  a  relationship 
for  the  critical  breakup  velocity  should  be  pursued. 

3.  The  ice  cover  utilized  in  this  study  served  only  as  a 
solid  boundary  condition.  The  effects  of  a  real  ice  cover 
should  be  explored. 
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TABLE  3.1 
OIL  PROPERTIES 


OIL  I  Rapeseed  oil 


Vi scosity 
Surface  Tension 
Density 


=  67.33  centipoise  at  21.1°  C 

=  14.71  dyne/cm 

=  0.915  gm/cm^ 


011  II  CWO  Base 


Vi scosity 

=  22.75 

centipoi se 

Surface  Tension 

=  28.24 

dyne/cm 

Density 

=  0.862 

gm/cm^ 

Heavy  Oil 

Vi scosity 

=  3794. 

centipoise  at  21.1°  C 

Surface  Tension 

=  32.33 

dyne/cm 

Density 

=  0.910 

gm/cm^ 

All  oils  were  properties  determined  prior  to  use  in  the  flume. 
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TABLE  5.1 
TEST  SCHEDULE 


RUN 

NUMBER 

ICE 

OIL 

VOLUME 

1  iter 

VELOCITY 
cm.  /s 

WATER 

DEPTH 

cm. 

BOOM 

DRAFT 

cm. 

LENGTH 

cm. 

Series  A 

1 

no 

I 

14 

24.1 

44 

9 

104 

2 

no 

I 

28 

24.1 

44 

9 

168 

3 

no 

I 

28 

25.5 

44 

9 

152 

4 

no 

I 

28 

27.0 

44 

10 

145 

5 

no 

I 

28 

28.3 

44 

10 

137 

6 

no 

I 

28 

29.7 

44 

12 

130 

7 

no 

I 

28 

27.6 

47 

10 

154 

8 

no 

I 

42 

24.1 

44 

9 

236 

9 

no 

I 

42 

25.5 

44 

10 

221 

10 

no 

I 

42 

27.0 

44 

12 

198 

11 

yes 

I 

42 

27.0 

44 

15 

191 

12 

yes 

I 

42 

25.5 

44 

15 

213 

13 

yes 

I 

42 

23.6 

46 

15 

251 

14 

no 

I 

42 

23.6 

46 

12 

267 

15 

no 

I 

28 

27.8 

47 

12 

152 

Series  B 

16 

no 

I 

21, 

24.1 

44 

10 

150 

17 

no 

II 

21 

22.2 

47 

8 

235 

18 

yes 

II 

21 

22.0 

47 

8 

213 

19 

no 

II 

21 

28.6 

44 

9 

190 

20 

yes 

II 

21 

28.6 

44 

9 

167 

21 

no 

II 

21 

26.1 

44 

10 

218 

22 

no 

II 

21 

24.1 

44 

10 

187 

23 

yes 

II 

21 

28.3 

44 

10 

169 

24 

yes 

II 

21 

28.3 

44 

10 

152 

• 
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TABLE  5.1  con't 


RUN 

NUMBER 

ICE 

OIL 

VOLUME 

liter 

VELOCITY 
cm.  /s 

WATER 

DEPTH 

cm. 

BOOM 

DRAFT 

cm. 

LENGTH 

cm. 

Series  C 

25 

no 

II 

21 

25.3 

44 

9 

213 

26 

no 

II 

21 

26.8 

44 

9 

179 

27 

no 

II 

28 

28.2 

44 

9 

213 

28 

no 

II 

28 

29.2 

44 

9 

192 

29 

no 

II 

28 

30.6 

44 

9 

178 

30 

no 

II 

28 

26.2 

44 

9 

229 

31 

yes 

II 

28 

26.2 

44 

9 

206 

32 

yes 

II 

28 

27.5 

44 

9 

190 

33 

no 

II 

28 

27.5 

44 

9 

210 

25 

no 

II 

21 

25.3 

44 

9 

213 

26 

no 

II 

21 

26.8 

44 

9 

179 

27 

no 

II 

28 

28.2 

44 

9 

213' 

28 

no 

II 

28 

29.2 

44 

9 

192 

29 

no 

II 

28 

30.6 

44 

9 

178 

30 

no 

II 

28 

26.2 

44 

9 

229 

31 

yes 

II 

28 

26.2 

44 

9 

206 

32 

yes 

II 

28 

27.5 

44 

9 

190 

33 

no 

II 

28 

27.5 

44 

9 

210 

34 

no 

II 

28 

28.5 

44 

9 

201 

35 

yes 

II 

28 

28.5 

43 

9 

182 

36 

no 

II 

28 

30.4 

43 

9 

165 

37 

no 

I 

21 

24.9 

45 

9 

137 

38 

yes 

I 

21 

24.9 

45 

9 

122 

39 

yes 

I 

21 

26.1 

45 

9 

114 

40 

no 

I 

21 

26.1 

45 

9 

122 

41 

no 

III 

14 

23.4 

45 

9 

91 

42 

yes 

III 

14 

23.4 

45 

9 

84 

43 

yes 

III 

14 

22.0 

45 

9 

91 

44 

no 

III 

14 

22.0 

45 

9 

99 

r 
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TABLE  5.2 

INTERFACIAL  VELOCITY  DATA 


RUN 

NUMBER 

ICE 

Uo 

cm 

s 

ui 

cm 

s 

h 

cm 

ai 

dyne 

2 

cm 

ci 

x  104 

X 

cm 

Rx 

x  105 

16 

no 

24.1 

12.0 

4.5 

1.76 

6 

83 

2.00 

17 

no 

22.2 

9.0 

2.5 

2.38 

9 

160 

3.55 

18 

yes 

22.0 

10.0 

3.0 

2.93 

12 

138 

3.04 

19 

no 

27.3 

18.0 

2.9 

4.10 

11 

115 

3.14 

20 

yes 

27.5 

17.0 

3.4 

4.40 

12 

92 

2.51 

21 

no 

22.5 

16.0 

2.6 

4.06 

13 

138 

3.52 

21 

no 

25.5 

16.0 

2.4 

4.40 

13 

83 

2.12 

22 

no 

28.3 

15.5 

3.0 

3.41 

8 

132 

3.74 

23 

yes 

28.3 

14.0 

3.8 

3.23 

8 

114 

3.22 

24 

yes 

29.3 

20.0 

3.7 

4.75 

11 

77 

2.26 
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►  U  <  0.2  m/s 

a)  LOW  VELOCITY 


near  boom 


b)  MEDIUM  VELOCITY 


c)  HIGH  VELOCITY 


FIGURE  1.1  CHARACTERISTIC  OIL  SLICK  PROFILES 
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a)  Headwave  Definition 


b)  Viscous  Zone  Definition 


FIGURE  2.1 


DEFINITION  SKETCH 
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FIGURE  2.2 


VISCOUS  ZONE  CONTROL  ELEMENT 
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a)  Forces  acting  on  a  differential  element  of  length  of  the  oil 
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FIGURE  2.3 


MILGRAM  AND  VAN  HOULTEN  OIL  SLICK 
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FIGURE  2.4  INSTABILITY  DEFINITION  SKETCH 
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FIGURE  3.1  SCHEMATIC  FLUME  CONFIGURATION 
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VELOCITY 


- ► 


a)  ICE  FREE  CASE 


VELOCITY 


-  - *■ 


b)  ICE  COVERED  CASE 


FIGURE  4.1 


IDEALIZED  VELOCITY  DISTRIBUTION 
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A  X 


p0  g(h+Ah) 

P0  g(m  +  Am) 


a)  ICE  FREE  FORCE  BALANCE 


AX 
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Pw  9 
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p„g  (m  +  Am) 


m  =  h  +  q  ,  Am  =  Ah 


b)  ICE  COVERED  FORCE  BALANCE 


FIGURE  4.2  FORCE  BALANCE  CONTROL  ELEMENT 
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v 


O  Rapeseed  oil  ice  free 
•  Rapeseed  oil  ice  covered 
□  CWO  base  ice  free 
■  CWO  base  ice  covered 
v  Heavy  oil  Ice  free 
t  Heavy  oil  ice  covered 
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FIGURE  5.1 


EFFECTIVE  LENGTH 
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FIGURE  5.2  ICE  FREE  VELOCITY  DISTRIBUTION 
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FIGURE  5.3  ICE  COVERED  VELOCITY  DISTRIBUTION 
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FIGURE  5.4  CORRELATION  OF  INTERFACIAL  VELOCITY  WITH  FREE 
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FIGURE  5.5  CORRELATION  OF  SHEAR  STRESS  WITH  INTERFACIAL  VELOCITY 
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FIGURE  5.6  CORRELATION  OF  FRICTION  COEFFICIENT  WITH  REYNOLD  NUMBER 
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OIL  SLICK  PROFILES 
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